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ABSTRACT 

We have identified and studied a sample of 151 FR lis found in brightest cluster galaxies (BCGs) in the 
MaxBCG cluster catalog with data from FIRST and NVSS. We have compared the radio luminosities and pro- 
jected lengths of these FR lis to the projected length distribution of a range of mock catalogs generated by an 
FR II model and estimate the FR II lifetime to be 1.9 x 10 8 yr. The uncertainty in the lifetime calculation is 
a factor of two, due primarily to uncertainties in the ICM density and the FR II axial ratio. We furthermore 
measure the jet power distribution of FR lis in BCGs and find that it is well described by a log-normal distri- 
bution with a median power of 1.1 x 10 37 W and a coefficient of variation of 2.2. These jet powers are nearly 
linearly related to the observed luminosities, and this relation is steeper than many other estimates, although it 
is dependent on the jet model. We investigate correlations between FR II and cluster properties and find that 
galaxy luminosity is correlated with jet power. This implies that jet power is also correlated with black hole 
mass, as the stellar luminosity of a BCG should be a good proxy for its spheroid mass and therefore the black 
hole mass. Jet power, however, is not correlated with cluster richness, nor is FR II lifetime strongly correlated 
with any cluster properties. We calculate the enthalpy of the lobes to examine the impact of the FR lis on the 
ICM and find that heating due to adiabatic expansion is too small to offset radiative cooling by a factor of at 
least six. In contrast, the jet power is approximately an order of magnitude larger than required to counteract 
cooling. We conclude that if feedback from FR lis offsets cooling of the ICM, then heating must be primarily 
due to another mechanism associated with FR II expansion. 

Subject headings: cooling flows — galaxies: active — galaxies: clusters: general — galaxies: evolution — 



galaxies: jets 



1. INTRODUCTION 



It has long been known that fewer extremely high-mass 
galaxie s are observed than s imple, theoretical calculations 
predict (White & Frenkll991l) . As random overdensities from 
the big bang collapse, hot gas in the largest overdensities 
cools and condenses toward their centers to form the progen- 
itors of brightest cluster galaxies (BCGs). In the absence of 
any source of heating, gas from the protocluster will continue 
to radiatively cool onto the BCG at rates as high as ^ 1000 
M Q yr" 1 in so-called "cooling flows" (see lFabianlfl994l for a 
review). However, X-ray observations indicate that the gas 
infall rate at the centers of most clust ers is about 10% of 
that predicted by c ooling flow models dPeterson et al.ll200H 
iTamuraet al]|2001l) . 

It is clear that some mechani sm heats gas on the cluster 
scale roughly isotropically (see iMcNamara & Nulsenl 120071 
for a recent review). The best candidate is some form of radio- 
mode feedback. Since BCGs tend to be disproportionately 
radio-loud, it has been suggested that most BC Gs go through a 
radio-loud phase dBurn si ll 990; Best et al. 2005). In this phase 
an AGN emits two radio jets which inflate large (^100 - 1000 
kpc), overpressured lobes in the intracluster medium (ICM). 
These radio lobes are often spatially coincident with X-ray 
cavities, implying that they physically displace a nd inject sub- 
stantial energy into some reg ions of the ICM dFabian et al.l 
2000; McNamaraetaDIIOO^. The power required to inflate 
these cavities appears to be sufficient to counteract c ooling 
of the hot gas dBirzan et al.ll2004t iRafferty et al"1l2006l) , Nev- 
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ertheless, it remains unclear how this energy heats the ICM 
sufficiently isotropically. Shocks and sound waves associated 
with the expansion of the lob e can transfer some heat from 
the lobe to the general IC M dJones et al.ll2002t iForman et al.l 
2005 UFabianetaill2005h a s could cosmic ray diffusion from 
the lobes (Mathews 2009). Conduction can ameliorate this 
problem to some extent, but to properly estimate the heating 
due to cond uction requires d etailed magnetohydrodynamic 
simulations dDolag et alj|2004l) . 

Large, extended radio sources have been well-studied 
for many decades. An early morphological classification 
scheme of doub l e-lobe d radio sources was developed by 
Fanaroff & Rilevl (119741) . Core-brightened sources are clas- 
sified as Type I sources (FR Is) and edge-brightened sources 
are classified as Type II sources (FR lis). Although the clas- 
sification is purely morphological, FR lis are generally more 
luminous than FR Is. The jet power supplied by the AGN is 
believed to play an important role in the FR I/II dichotomy, 
but the occasional cases of "hybrid" sources which exhibit 
one FR I lobe and one FR II l obe indicate that environmenta l 
factors are not negligible (e.g., Gopal-Krishna & Wiita 2000). 

Because FR lis have sharp edges, many important prop- 
erties (e.g., projected length, axial ratio, and luminosity) are 
well defined and independent of the sensitivity of the ob- 
servation. By contrast, FR Is exhibit long plumes of fad- 
ing radio emission, so properties like the length and total 
luminosity are more dependent on the sensitivity of the ob- 
servation. FR lis are therefore easier to study as probes 
of AGN activity and their effects on the ICM with radio 
observations alone. Furthermore, the advance speeds and 
axial ratios of FR lis are similar across several orders of 
magnitude in lobe length, indicating that FR II evolution 
is approximately self-similar. This inference has addition- 
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ally been supported by more detail ed hydrodynamical sim- 
ulations (Carva lho & O'Deal l2002allbh . The assumption of 
self-similarity has allowed the possibility of analytical or 
semi-analytical models to describe FR II evolution. A num- 
ber of such models have been developed over the past two 
decades. Among the most complete and most widely used 
is the model presented in iKaiser & Alexander! (1 19971) and 
iKaiser. Dennett- Thorpe. & Alexanderl d 1997b . henceforth re- 
ferred to as the KDA model. We use the KDA model exclu- 
sively throughout this paper. 

The calculation of the length evolution of the lobe is 
straightforward and is similar in all FR II models. Ear- 
lier FR II studies have used FR II lobe lengt h distributions 
to estimate FR II age s of 10 7 to 10 s years dBlundell et al.l 
119991: iBird etaill2008l) . These measurements are in reason- 
able agreement with spectral aging estimates and measure- 
ments of the buoyant rising times of ghost ca vities, which find 
FR II ages in the range of 10 6 to 10 8 yr dAllen et al .1120061: 
lO'Dea et alJl2009h . These estimates generally only reflect the 
current age of individual sources, or subsets of the population 
that may not be representative, and do not provide a direct es- 
timate of the lifetime or energy of typical sources. While it is 
possible to infer an FR II lifetime from an age distribution, it 
requires good knowledge of the selection efficiency as a func- 
tion of age. Studies which attempt to infer the FR II lifetime 
from length distributions can con s train their selec tion biases 
more easily dBlundell et all 1 19991: IBird et alJl2008h . but still 
suffer from uncertainties in the jet power. 

The FR II jet power is a difficult quantity to determine 
and the few FR lis for which it has been estimated have 
tended to be unusually bright, and hence, high-powe r sources 
dMachalski et alj 120041: iMcNamara & Nulsenl l2007h . Infer- 
ring the general FR II jet power distribution from the high- 
power end alone, however, requires substantial and uncer- 
tain extrapolation. Nevertheless, recent studies have begun 
to probe the low-power end of the FR II jet power distribution 
and have found typic al jet powers of 10 37 W and jet pow- 
ers as low as 10 3 ° W dCavagnolo et alJ bOlOt llto et al.ll2008l: 
iPunslv & Zhangl [201 lb . These estimates place the median 
FR II power lower than older estimates by at least an order 
of magnitude. Although FR II lifetime estimates do not scale 
sensitively with FR II power, a substantial overestimate of the 
jet power will lead to an underestimate of the age or lifetime. 

In this work we measure the ages and powers of 151 
FR lis in BCGs from the MaxBCG galaxy cluster catalog of 
iKoester et alj d2007al) . We simultaneously fit the length and 
luminosity of the FR lis with the KDA model and derive an 
estimate of the power distribution. We then obtain an im- 
proved lifetime and duty cycle estimate and examine the im- 
portance of FR lis in ICM heating. The paper is organized as 
follows: in ^2]we describe the sample selection algorithm, the 
luminosity and projected length measurements, and list gen- 
eral properties of the FR lis in our sample; in §|3]we motivate 
our use of the KDA model, provide a summary of the model's 
important features, and describe our technique to estimate the 
jet powers; in $4] we fit the observed length distribution to 
a set of mock catalogs generated by a Monte Carlo simula- 
tion of FR II observations; in $5] we investigate correlations 
between FR II power and lifetime with various cluster prop- 
erties and examine the FR II fraction as a function of cluster 
properties; in §|6]we compare our results with those in the lit- 
erature and discuss the heating effects of FR lis on the ICM; 
we summarize our results in §|7] 

Throughout the paper we assume a ACDM cosmology with 



n„, = 0.27, Q. A = 0.73, and H = 70 km s" 1 Mpc" 1 . 

2. SAMPLE SELECTION 

Our sample selec tion process is similar to the process used 
in lBird et al.l(l2008l) . namely we identify FR lis hosted by cen- 
tral galaxies in dense environments. The major difference is 
that we use the MaxBCG catalog of IKoester et alJ <2007al) to 
identify B CGs rather than the catalog of galaxy groups pub- 
lished bv lBerlind et al.l (12006b . We describe properties of the 
galaxy cluster catal og i n $2. II and properties of the catalogs 
of radio sources in $2.21 Our process for cross-correlating the 
two catalogs is described in 32.31 32.41 describes the sample 
selection method and properties of the final FR II sample. 

2.1. Cluster Catalog 

The MaxBCG catalog consists of 13,823 galaxy clusters se- 
lected from the Sloan Digital Sky Survey with the maxBCG 
red-sequence method. The red sequence method attempts 
to find BCGs by selecting galaxies lying on the E/S0 ridge- 
line (a region of color-magnitude space where old, passively 
evolving galaxies are found) which are also brighter than 
their neighboring galaxies. The algorithm identifies potential 
BCGs in an overdensity of galaxies with similar colors and 
redshifts. Th e red se quence method is described in detail in 
IKoester et al.l d2007al). T ests using mock catalogs performed 
by lKoester et al.l (l2007bh indicate that the MaxBCG catalog is 
at least 90% pure and 85% complete for clusters larger than 
10 14 M Q . 

The number of galaxies with similar colors and redshifts in 
the overdense region, N ga j, is an estimate of the total number 
of galaxies in the cluster. N ga \ is later refined to A^ 00 , the 
number of galaxies within the radius at which the galaxy den- 
sity is 200 times larger than the mean density of galaxies with 
-24 < M r < -16 mag. Throughout this paper we use A^ 200 as 
a proxy for the cluster richness. 

The MaxBCG cluster catalog spans a range in photomet- 
ric redshift of z = 0. 1 - 0.3 and a range in richness of A^ 200 = 
10- 188. 39% of the BCGs in the MaxBCG catalog also have 
spectroscopic redshifts and the dispersion between the spec- 
troscopic and photometric redshifts is 0.01. Throughout this 
paper we use spectroscopic redshifts whenever they exist and 
otherwise use photometric redshifts. The distribution of clus- 
ter richness and redshifts in the MaxBCG catalog is shown in 
Figure Q] 

Although the MaxBCG catalog does not provide cluster 
mass estimates, the clusters span a range of r-band luminosity 
from 7 x 10 10 L Q to 3 x 10 12 L Q with a median luminosity of 
2 x 10 11 LfT> For a typical clus t er mass- to-light rat io of ^350 
M / L dCarlberg et al.l[l997t ISheldon et al.ll2009h . the typ- 
ical cluster mass in the MaxBCG catalog is ^7 xlO 14 M Q 
with a spread of about one order of magnitude in either direc- 
tion. 

2.2. Radio Catalogs 

We selected radio sources associated with the BCGs of 
these galaxy clusters using the July 16, 2008 version of the 
FIRST radio catalog dWhite et al.lll997h . The FIRST radio 
catalog consists of 816,331 radio sources identified from 9055 
deg 2 observed as p art of the FIRST radio survey at 1 .4 GHz 
(Beck eret al.lll995h . Data for this version of the catalog were 
obtained with the Karl G. Jansky Very Large Array (VLA) 
from 1993 through 2004. The angular resolution of the FIRST 
survey is 5. "4 and the typical rms flux noise is 0.15 mJy. Dirty 
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FIG. 1 . — The distribution of cluster richness and redshift in the MaxBCG catalog (small gray points), the full FR II sample (large open and filled circles), and 
the volume-limited FR II sample (filled circles alone). Histograms of the cluster richness and redshift are also shown. The unfilled area represents the counts for 
all clusters in the MaxBCG catalog, the light gray area the counts for those clusters with FR lis in our full sample, and the dark gray area the counts for those 
clusters with FR lis in our volume-limited sample. 



beams from the VLA were CLEANed to produce images for 
the FIRST survey with 1 ."4 pixels. 

The FIRST radio catalog was produced from the survey 
data using a specialized source extraction algorithm dubbed 
HAPPY. HAPPY searches for local maxima and attempts to 
fit them to up to four elliptical Gaussian components. Fits 
which produce physical results and are not too close to an im- 
age edge are accepted into the catalog. The catalog is then 
refined to eliminate duplicate and spurious sources. To be 
considered a detection, a source must be at least five times 
larger than the rms noise. The minimum flux of any source is 
therefore typically above 0.75 mJy. 

The resulting catalog provides, among other parameters, 
the coordinates, the local rms noise, the peak and integrated 
flux, and a warning flag if the source is a potential sidelobe 
of a nearby bright source. The two flux measurements are 
correc ted for CLEAN bias dBecker et al.lll995l ICondon et ail 
[19981) . 

Because the FIRST survey has relatively high resolu- 



tion and treats extended sources as collections of point-like 
sources, much of the flux of extended objects can be resolved 
out and missed, thereby leading to underestimates of the to- 
tal flux. Since FR lis are extended, we complement our flux 
measurements of sources from the FIRST catalog using the 
NRAO VLA Sky Survey (NVSS). NVSS is a similar radio 
survey to FIRST, but has an angular resolution of 45" rather 
than 5 "4, thereby pro viding more accurat e flux measurements 
for extended sources (ICondon e t al. 1998). The flux limit for 
the NVSS catalog is S ~ 2.5 mJy. Although this is larger than 
FIRST by a factor of approximately three, extended sources 
are broken into fewer components in the NVSS catalog and 
we lose no FR lis in the NVSS catalog. 

2.3. Cross-Correlating the Two Catalogs 

We first selected a sample of FR II candidates by identifying 
radio sources between 10 "8 and 3.4' of every BCG The upper 
bound of 3.4' was chosen to be small enough to limit the num- 
ber of spurious candidates, yet large enough that real FR lis 
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would not be eliminated. At a redshift of z = 0.1, the low- 
est in our sample, this cut would eliminate any FR lis longer 
than 370 kpc. Since we found no FR lis longer than 287 kpc, 
we conclude that this choice of maximum angular separation 
has no effect on our sample size. The lower bound of 10."8 
was chosen so that the FR lis would be at least a factor of 
two greater than the FIRST resolution. The selection function 
for short FR lis near the FIRST resolution limit is difficult 
to characterize; some FR lis near the FIRST resolution limit 
were discovered by our selection algorithm, whereas others, 
due to the chance peculiarities of their morphology or sur- 
roundings, may have been missed. By selecting only FR lis 
whose lobes each subtend at least two resolution elements of 
FIRST, we ensure that our selection algorithm only identifies 
FR lis that are resolved by FIRST. 

Except in rare cases when a radio jet is oriented very nearly 
along the line of sight and has one lobe strongly beamed to- 
wards the Earth and the other lobe strongly b eamed away, ra - 
dio jets appear to be extremely symmetrical (Sch euerlll995l) . 
We thus identified FR II candidates as radio sources that were 
equidistant from the central BCG to within 30% and at least 
135° away from each other relative to the BCG. These selec- 
tion criteria yielded an initial sample of 617 FR II candidates. 
We then superimposed FIRST contours on SDSS images of 
these FR II candidates and inspected them by eye to eliminate 
contaminants. There were two major sources of contamina- 
tion: (1) pairs of radio galaxies on opposite sides of the BCG 
or (2) a radio galaxy paired with a faint radio source unasso- 
ciated with any optical counterpart. Aside from their spatial 
coincidence with optical sources, contaminants in both cases 
could be easily distinguished from faint FR II lobes due to the 
relatively large asymmetry between the position of the two 
sources with respect to the BCG, the relatively large asymme- 
try in flux, and the lack of any extended emission toward the 
BCG. We also removed 31 obvious FR Is from the sample. 
This resulted in a preliminary sample of 432 FR lis. Less- 
obvious FR I sources were removed in a more quantitative 
second pass described in the following section. 

2.4. FR II Classification and Sample Selection 

Although FR Is can often be distinguished from FR lis 
based on luminosity alone, the formal divisi on between FR I 
and F R II sources is purely morphological (Fanaroff & Rilevl 
fl97l . Sources for which the hotspot is closer to the tip of 
the lobe are classified as FR lis, and sources for which the 
hotspot is closer to the central galaxy are classified as FR Is. 
Except in extreme cases, since projection effects cannot dras- 
tically change the ratio between the distance to the hotspot 
and the distance to the lobe tip, the projected positions of the 
hotspot, lobe tip, and central galaxy are sufficient to eliminate 
any remaining FR Is from our sample. 

We therefore measured the position of the hotspot and the 
tip of the lobe for each FR II in the sample. We identified the 
position of the hotspot as the centroid at the local maximum 
near the FR EE's edge. It is inherently more difficult to deter- 
mine the tip of the lobe. Although the lobes of FR II sources 
have better-defined edges relative to FR Is, the tip of a lobe 
is nevertheless not entirely well-defined due to variations in 
the sensitivity of the radio data and the often complex mor- 
phology of the source. We defined the tip of the lobe as the 
most distant pixel from the central BCG that was contiguous 
with the main lobe and above the FIRST detection threshold 
by five times the rms noise. Once we determined the position 
of the hotspot and tip of the lobe, we removed less-evident 



FR Is from the sample. 

We derive the projected physical length of each source from 
the separation of each lobe tip from the BCG and the angular 
diameter distance to the cluster. The errors on the projected 
source sizes are dominated by random errors in the photomet- 
ric redshift and the uncertainty in the definition of the lobe 
tip. We estimate that the total uncertainty in the physical lobe 
length is < 5%. As shall be shown in § 34.2l and l4.3l the uncer- 
tainty in our lifetime calculation is dominated by uncertainty 
in the ICM and the FR II model, rather than the uncertainty in 
physical lobe length. 

We obtained the fluxes for the FR lis from the NVSS cata- 
log. Because NVSS, like FIRST, breaks up extended sources 
into discrete components, the flux from each component was 
summed to give the total flux from each side of the FR II 
source. Since NVSS cannot resolve the double-lobed mor- 
phology of the smaller sources in the sample, it was often 
impossible to assign a separate flux to each lobe. In these in- 
stances we measured the total flux from both lobes. When 
fluxes from individual lobes were needed, such as when esti- 
mating the power and age of the jet as described in ^3.31 we 
simply assigned each lobe half the total flux. Since this pre- 
scription requires that both lobes be FR lis, we made a final 
cut in our sample to eliminate any radio sources that consisted 
of an FR II lobe paired with an FR I lobe. This cut eliminated 
42 nominal FR I/II mixed sources. We note that the vast ma- 
jority of these mixed sources are not "true" FR I/II hybrid 
sources in the usual sense of a well-defined FR I plume on 
one side and a well-d efined FR II lobe on the other (so-called 
"HYMOR" sources. iGopal-Krishna & Wlt3l2000h . Rather, 
in these sources the hotspots of both lobes are approximately 
halfway along the length of the lobe and one hotspot is slightly 
beyond halfway to the lobe tip and the other slightly before. 
Application of all of these selection criteria resulted in a flux- 
limited sample of 151 FR II lobe pairs. We perform the bulk 
of analysis with this sample and will simply refer to it as "the 
sample" or the "full sample" rather than the flux-limited sam- 
ple. The distribution of the FR lis in our sample in redshift 
and host cluster richness is shown in Figure Q] and the pro- 
jected length distribution is shown in Figure|2] 

This sample of 151 FR II lobe pairs suffers from several 
significant selection biases. Specifically, intrinsically longer 
jets are easier to resolve at larger distances, and more lumi- 
nous FR lis are easier to detect at larger distances. Since 
more powerful FR lis are both longer and more lumino us at a 
given age than less powerful FR lis dKaiser et alJI 19971) . this 
sample is biased towards more powerful FR lis at higher red- 
shifts. These biases become problematic when attempting to 
determine whether FR II lifetime or duty cycle change with 
redshift. To ameliorate this problem, we introduced a cut into 
our sample that required every FR II source meet our selec- 
tion criteria at z = 0.3, the largest redshift of the clusters in 
our sample. That is, if the source were at z = 0.3, the angu- 
lar separation between the two lobes of the jet would have 
to be at least 10."8, and the flux from the FIRST component 
for both hotspots would have to be above the FIRST detec- 
tion limit of five times the rms background noise. There are 
108 pairs of FR II lobes in this sample, which we will refer to 
as the volume-limited sample. The length distribution of the 
volume-limited sample is also shown in Figure [2] 

There is likely a slight overcorrection in the volume-limited 
sample, but one that is difficult to quantify. Because ex- 
tended sources appear more compact at larger distances, high- 
resolution surveys like FIRST can collect more flux from the 
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FIG. 2. — The length distribution of the full (solid line) and the volume- 
limited (dashed line) FR II samples. Since shorter sources are not resolved 
by FIRST at high redshifts, the volume-limited sample contains fewer short 
FR lis than the flux-limited sample. 

source. Thus if an extended source at z < 0.3 is moved out to 
Z = 0.3, the flux detected by FIRST would be slightly larger 
than one would expect based on the change in distance alone. 
Since this effect is difficult to characterize accurately due to 
the complex morphology of the sources, we simply note that 
low-redshift sources in our volume-limited sample are slightly 
biased towards larger lengths. 

The properties of all of the FR II sources and their host 
clusters are listed in Table Q] Figure [3] displays examples of 
several randomly selected FR lis from our sample and Figure 
|4]presents a P-D diagram of the full sample. The sources span 
a range of lengths of 48 - 287 kpc with a median length of 95 
kpc. The host clusters span a range of richness from 10 (the 
minimum richness in the MaxBCG catalog) to 69 with a me- 
dian richness of 14 and the full redshift range of the MaxBCG 
catalog of z = 0.1-0.3. The sample spans a range in specific 
luminosity of 4.4 x 10 23 -7.2 x 10 25 W m" 2 Hz" 1 with a me- 
dian luminosity of 8.0 x 10 24 W m" 2 Hz" 1 (all at an observed 
frequency of 1.4 GHz). Note that we do not ^-correct the lu- 
minosities of the FR II sample at this point, but acco unt for 
A"-corrections in the mock catalogs (described in 33. 5t . 

These luminosities are smaller by at least an order of mag- 
nitude than th e luminosities of FR II sources in many earlier 
studies (e.g.. lLaing et alj 119831; ISubrahma nvan et al.1 119961; 
ICotter et al.l ll996) and are smaller than the luminosities of 
well-studied FR lis like Cygnu s A, 3C 47, and 3 C 295 by 
over three orders of magnitude dBraude et al.lll969l) . Indeed, 
75% of the sources in ou r sample fall below the FR I/II de- 
marcation line set out in lLedlow & ()\venl(ll996). While it is 
possible that some sources close to our resolution cutoff may 
be FR Is confused for FR lis, the vast majority of our sources 
are well resolved. Specifically, ^75% of our sample has an 
angular separation of at least three FIRST beams, ~50% of 
our sample has an angular separation of at least four FIRST 
beams, and ~25% of our sample has an angular separation of 
at least six FIRST beams. 

The relatively low luminosity of the FR lis in our sample is 
due to the FIRST survey's superior sensitivity and resolution 
with respect to earlier radio surveys. Probing the faint end of 



the FR II luminosity function allows us to better constrain the 
lower-limit of FR II jet power and study the characteristi cs of 
more typical FR lis. We address this subject in detail in 33.31 

3. MODELING THE LENGTH DISTRIBUTION 

We estimate the lifetime and jet powers of FR lis with de- 
tailed models of their luminosity evolution, environments, and 
the sample selection. This approach requires a good model of 
FR II leng th an d luminosity evolution, and our choice is dis- 
cussed in 33.11 The model for the ICM can impact the lifetime 
determination as well, so we discuss the properties of the ICM 
models we use in 33.21 In 33.31 we introduce the method by 
which we estimate FR II jet powers, a major contribution to 
the uncertainty in previous FR II lifetime measurements. We 
present a relationship between the jet power and the radio lu- 
minosity in 33.41 We finally discus s ou r use of mock catalogs 
to account for selection effects in 33.51 

3.1. Properties of the FR II Model 

To determine the FR II source lifetime we need an accu- 
rate, analytic model of the length and luminosity evolution of 
the sources. A number of analytic mo dels of FR II source 
evolution exist in the literature (e.g.. [Bicknell et al.1 fl997l; 
Kaiser et al.l 119971; iBlundell et alj [l999t iManolak ou & KirU 
2002t iKino & Kawakatull2005h . The Bicknell, KDA, Blun- 
dell, and Manolakou models all assume self-similar growth of 
the jet. Although this assumption appears to be approximately 
correct, hydrodynamical simulations by ICarvalho & O'Deal 
( 2002a, b) show that it fails in detail in low Mach number jets 
in ICMs with r elatively flat dens i ty pro files. The more re- 
cent model by IKino & Kawak atu ( 2005) avoids the assump- 
tion of self-similarity, although it only describes the total ki- 
netic power rather than the radio luminosity of the source. 
Since the density profile of the ICM of galaxy clusters tends 
to be steep and the jets in FR lis have large Mach numbers, 
the self-similarity assumption is not problematic for our pur- 
poses. The radio luminosity of the FR II sources plays an im- 
portant role in our determination of their lifeti mes, so we find 
the Ki no model to be unsuitable for our study. Barai & Wiital 
(2006) compared the Bicknell, Manolakou, and KDA models 
to FR II radio galaxies in the 3CRR, 6CE, and 7CRS radio 
surveys and found that the KDA model best matched data in 
the P-D-z-a plane. Based on this r esult, and for pu rposes of 
comparison with a similar study by Bir d et all (120081) . we use 
the KDA model to describe the length and luminosity evolu- 
tion of jets in our sample. 

The KDA mode l was re cently reworked into a simpler form 
by iKaiser & Bestl (|2007); we adopt their notation for this 
paper. We note, however, that we have made one modi- 
ficat ion to the original KDA m odel. In a subsequent pa- 
per, |KaiieE&Alexandei 09992) performed hydrodynamical 
simulations to test the accuracy of the KDA model. While 
the analytic model agreed with hydrodynamical simulations, 
Kaiser & Alexander foun d that an approximation made in 
IKaiser & Alexander! d 19971) of 

^=A 2 (1) 

Pc 

overestimates /?/,/ p c , where pf, is the pressure at the head of 
the jet, p c is the pressure within the cocoon, and A is the axial 
ratio. Kaiser (2000) provided an empirical parameterization 
of this ratio based on earlier simulations, giving 

/ . 2.04-0.25/3 

^ =(2.14-0.52/3) (-J , (2) 
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FIG. 3. — Twelve randomly chosen FR II radio sources from the sample. These are SDSS r-band images with contours from the FIRST survey superimposed. 



TABLE 1 
FR II Sample 



SDSS ID 


z 


N200 


BCG 


L z , BCG 


Fnvss 




h 


h 


h 


h 


In volume- 








(1O IO L ) 


(1O 1O L ) 


(mJy) 


(10 24 WHz- 1 ) 


(kpc) 


(kpc) 


(arcsec) 


(arcsec) 


limited sample? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


J001247.6+004715.8 


0.154 


14 


5.46 


6.83 


61.7 


3.9 


46.9 


56.0 


13.2 


15.8 


No 


J003614.2-090451.8 


0.270 


10 


4.43 


5.24 


17.0 


3.8 


121.0 


127.8 


18.1 


19.1 


Yes 


J004312.9-103956.1 


0.128 


19 


7.71 


9.60 


148.4 


6.2 


119.6 


122.8 


41.4 


42.5 


Yes 


J013157.8-081955.0 


0.140 


11 


7.51 


9.65 


469.0 


24.2 


55.5 


85.9 


17.3 


26.8 


No 


J073826.2+451719.1 


0.221 


10 


6.36 


8.07 


91.9 


13.0 


65.3 


80.3 


12.3 


15.1 


No 


J074616.9+220203.5 


0.260 


19 


6.16 


7.51 


30.4 


6.2 


86.1 


110.1 


13.5 


17.2 


Yes 


J075457.8+210129.0 


0.256 


14 


5.42 


6.78 


21.7 


4.3 


64.2 


74.0 


10.8 


11.8 


Yes 


J075614.8+251340.4 


0.202 


11 


8.67 


10.77 


188.8 


21.7 


75.5 


111.1 


15.7 


23.2 


Yes 


J080107.0+175845.3 


0.146 


24 


5.05 


6.52 


890.7 


49.9 


78.1 


104.8 


23.3 


31.3 


Yes 


J080641.4+494628.4 


0.245 


13 


5.90 


7.66 


84.2 


15.0 


84.9 


112.4 


14.2 


18.8 


Yes 



NOTE. — FR II sources in in our sample ordered by right ascension. (1): SDSS identifier. (2): BCG redshift. (3): Cluster richness. (4, 5): The r- 
and z-band luminosities of the BCG, respectively. (6): Total flux as measured by NVSS in mJy. (7): Specific radio luminosity of the FR II calculated 
from the NVSS flux in 10 24 WHz _1 . (8, 9): Projected length of the two lobes in kpc. (10, 11): Angular separation of the two lobes in arcseconds. (12): 
Identifies whether the FR II is also in the volume-limited sample. Table[T]is presented in its entirety in the electronic edition. A portion is shown here 
for guidance regarding its form and content. 



where (3 is the power law slope of the ICM density profile. 
This modification has the following effect on the KDA model. 
The length of the lobe in the original KDA model is given by 



where C[ is 



C] 



D = c, 



Q_ 



1/(5-/3) 



,3/(5-/3) 



(r,+ ixrw)(5-/3) 3 



1 8tt 9[r, + (T, - 1 )A 2 /2] - 4 - /3 



1/(5-/3) 



(3) 



(4) 



t is the age of the FR II, Q is the jet power, and T x and T; are 
the adiabatic indices of the ICM and the lobe, respectively. 
By substituting Equation (Q]) for Equation (O, the constant c\ 
becomes 



J-/3 . 



(2.14-0.52^)A 2 



18tt 



2.04— 0.25/3 ^ 



(r t +i)(r,-i)(5-/3) 3 



(5) 



,9[r, + (r,-l)A 2 /2]-4-/3 
For our choice of f3 = 1.9 (discussed in 33.2b . this modifica 
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FIG. 4. — Radio luminosity vs. lobe length for FR lis in the full (open and 
filled circles) and volume-limited samples (filled circles only). There is no 
clean line in length-luminosity space that divides the two samples because 
the sample cut was made with the length to the hotspot, while the length to 
the edge of the lobe is shown in this figure and the ratio between the hotspot 
and lobe distance can vary by up to a factor of two. For comparison, evolu- 
tionary tracks from the KDA model for four extreme cases are shown. The 
solid lines show an ICM model with a = 260 kpc and p = 5.9 X 10~ 24 kg/m 3 , 
characteristic of a rich cluster, and the dashed lines show an ICM model with 
a = 2 kpc and p = 7.2 X 10~ 22 kg/m 3 , characteristic of a small group. Both 
the upper solid and dashed lines represent a jet with a power of Q = 10 37 W, 
and the lower solid and dashed lines represent a jet with a power of Q = 10 36 
W. The left and right ticks of the solid lines represent ages of 10 s and 3 X 10 s 
yr, respectively, and the left and right ticks of the dashed lines represent ages 
of 3 X 10 7 and 10 s yr, respectively. 

tion shortens the predicted lengths by approximately 40% and 
increases the age by approximately 70%. 

The salient feature of the KDA model is that it predicts 
the source length and specific luminosity of a jet given an 
age, jet power, redshift, and a number of additional param- 
eters that describe the jet's shape and environment. The full 
list of inputs and several of the most relevant outputs of the 
KDA model is displayed in Table[2] along with the range over 
which we vary the parameters and our default values in this 
paper. Because the KDA model must be run several million 
times for a given choice of inputs in order to accurately model 
the length distribution, it is too computationally expensive 
to sample the entire range of p arameter space that the KDA 
model offers. iBird et all ([2008 ) found, however, that most of 
the parameters have only a negligible impact on FR II life- 
time calculations over physically plausible ranges. We there- 
fore limit our exploration of parameter space to only those pa- 
rameters which Bird et al. (2008) found to have a substantial 
impact on lifetime calculations, namely the density parameter 
pa" and the axial ratio A. 

3.2. Treatment of the ICM 

Apart from an assumption of self-similar growth, the KDA 
model also assumes that the ambient medium is described by 
a power-law profile of the form 



where a is the scale length of the distribution and p is the 
density at the scale length. As every other equation in the 



KDA model involving p and a incorporates this one, these 
two parameters are not independent, but instead enter into the 
model only in their combined form, pa 13 (/3, however, enters 
the mo del in separate equa tions and so is independent from a 
and p). iKaiser & Bestl (120071) term the product pa 13 the "den- 
sity parameter" and we adopt that nomenclature in this paper. 

The KDA model uses this simplified profile to make the 
calculation of more complicated quantities feasible, although 
at the expense of some shortcomings in the model. For ex- 
ample, while the KDA model adopts a density distribution for 
the ambient medium, there is no model for the pressure or 
the temperature. As a result, the KDA model tacitly assumes 
that the jet cocoon remains overpressured throughout its en- 
tire lifetime. Indeed, the KDA model require s that the lobe 
drive a strong shock into the external medium. IKaiser & Bestl 
(120071) point out that the pressure in the cocoon of an FR II 
source decreases over time, and argue that if the cocoon pres- 
sure falls below the ambient pressure at any point along the 
surface of the cocoon, Kelvin-Helmholtz instabilities will mix 
cooler gas from the external medium into the cocoon. This 
process would disrupt the jet and transform it into an FR I. 
If the pressure profile can be described by a flat core with a 
power-law tail, then low-power jets will become underpres- 
sured and transform into FR Is before escaping from the core 
region; however, more powerful jets will escape the core re- 
gion and will remain overpressured (a nd hence FR lis) for 
their entire lifetimes. With this model, Kaiser & Best! (120071) 
present an order-of-magnitude calculation and find that the 
critical power is 3 x 10 37 W. 

Given the lengths of the FR lis in our sample, the KDA 
model predicts that even FR lis with this minimum power are 
much more luminous than the sources in our sample. As de- 
scribed next in 33.31 to reproduce jets with roughly the same 
lengths and luminosities as the jets in our sample, the KDA 
model requires that the typical jet power be over an order of 
magnitud e smaller than t he mi nimum FR II jet power power 
found by IKaiser & Bestl d2007h . This would imply that the 
jets are underpressured prior to leaving the flat core of the 
ICM. There are two problems with this interpretation, how- 
ever. The first is that a detailed examination of the pressure 
profiles of galaxy groups and clusters reveals that, while com- 
plex and varied in shape, they g enerally have cusps a t their 
centers rather than flat cores (e.g- IVikhlinin et al.ll2006l) . This 
suggests that many underpowered jets would start their lives 
as FR Is rather than later disrupting from FR lis. Because 
the ICM density profiles measured by Vikhl inin et al.l (120061) 
can vary substantially from one group or cluster to the next, 
there is likely no single, universal critical power. Furthermore, 
because the ICM profile can exhibit plateaus and bumps, it 
is difficult to predict whether an FR II will ever become un- 
derpressured, and if so, at what radius, without knowing the 
detailed ICM profile of a group or cluster. The second prob- 
lem with the interpretation of the jets in our sample as un- 
derpressured is that the KDA model calculates the pressure at 
the head of the jet and assumes the pressure in the cocoon to 
be uniform. Disruption of the jet is likely to occur near the 
base of the jet, however, where the ICM pressure is greatest. 
Since the jet expansion speed can be comparable to, or even 
larger than, the sound speed within the cocoon, the pressure 
within the cocoon at the head of the jet will not necessarily 
be a good proxy for the pressure within the cocoon near the 
center of the FR II source, particularly towards the end of the 
jet's life. Indeed, hydrodynamical simulations suggest that the 
pressure within the cocoon varies along and across the axis of 
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TABLE 2 

Input & Output Parameters of the KDA Model 



r arameter type 


Parameter 


Default value 


Minimum value Maximum value 


■ • 
Description 


Primary input 


t 






FR II age 


Primary input 


Q 






Jet power 


Primary input 


z 






Redshift 


Secondary input 


A 


4.0 


2.0 16.0 


Axial ratio (see 34.2 tfor our definition of A) 


Secondary input 


P 


7.2 X 10~ 26 kg irT 3 


7.2 x 10" 26 kg nT 3 7.2 x 10~ 22 kg m" 3 


Density at uq 


Secondary input 


Ll{) 


391 kpc 


2 kpc 391 kpc 


Scale radius of ICM density profile 


Secondary input 


/3 


1.9 




Power-law slope of ICM density profile 


Secondary input 


m 


2.14 




Power-law slope of injection energy of particles in iet 


Secondary input 


V 


1.4 GHz 


... 


Observation frequency 


Secondary input 


Tx 


5/3 




Adiabatic index of the ICM 


Secondary input 


r, 


4/3 




Adiabatic index of the lobe 


Secondary input 


r a 


4/3 




Adiabatic index of the magnetic field energy density 


Secondary input 


Tmin 


1 




Minimum Lorentz factor for electrons in the jet 


Secondary input 


""/max 


10'° 




Maximum Lorentz factor for electrons in the jet 


Secondary input 


k 







Ratio of energy stored in non-radiating particles 
to the sum of the energy in the magnetic field and the 
relativistic electrons. 


Output 


i 






FR II Lobe length 


Output 


L 






FR II radio luminosity 


Output 


Pc 






FR II lobe pressure 


NOTE. — We designate input parameters which are 


likely to be significantly different among FR lis to be "Primary inputs" and therefore assign no 



default values to them. A similar study by Bird et al. (2008) found that variation of most of the secondary inputs does not substantially affect the result of the 
lifetime calculation. We therefore only varied those parameters which have a substantial effect on the lifetime calculation. Maximum and minimum values 
for parameters that we did not vary are given by ellipses. 



an FR II jet and that the cocoon can become underpressured 
with respect to the ambien t medium without disrupting the jet 
dCarvalho & O'Deal [2002b). Based on these considerations 
and our observations , we conclude t hat th e minimum FR II 
jet power derived by iKaiser & Best! (120071) is too high by at 
least an order of magnitude. 

The density profile of the ICM in galaxy clusters is diffi- 
cult to measure and the core den sity and scale radius can vary 
by several orders of magnitude (IMohr et al.lll999t iMul chaev 
l2000t|Vikhlinin et al.l2006tlFreeland & Wilcotsl201 lb . ICotteri 
( 1996) and Kais er et all (11997) found that the power-law slope 
of the ICM density profiles in galaxy clusters hosting FR lis 
are typically consistent with a value of f3 = 1 .9. We adopt this 
slope throu ghout this p aper. Larger values of /3 cannot sup- 
port FR lis ( Falla [l99ll) and smaller values result in a smaller 
variation in the density parameter in different ICMs, thereby 
makin g our lifetime res ults less sensitive to the ICM prop- 
erties. Bird et al.l (120081) demonstrated in their appendix that 
the lifetimes of FR lis depend upon the cube root of the den- 
sity parameter, so a three order-of-magnitude variation in den- 
sity parameter manifests itself as an order-of-magnitude un- 
certainty in the lifetime calculation of the jet. Although this 
dependence is weak, we note that quantities which depend on 
the volume of the lobe (e.g., the energy contained in the lobe) 
depend linearly on the density parameter. 

Insufficient X-ray data exist to measure the ICM density 
profiles for these clusters. Instead we calculate the FR II life- 
time separately for five different ICM density estimates drawn 
from the literature and then use optical richness as a proxy for 
ICM density in $5] To facilitate reference to each of the den- 
sity estimates we label each with the symbol \i- These five 
density estimates are listed in full in Table|3]with references, 
but we describe them in order of increasing de nsity parameter 
here a s follows. \\ lS the den sity assumed b y Blunde ll et al] 
(119991) . which was taken from Garringt on & Conwayl (1991 ) 
and assumed to be typical of a poor group, xi lS the de- 



fault density assumed by Kaiser et al. ( 1997) and was taken 
from Canizares et al. (f!987ir who a lso found it to be t ypical 
for poor groups. Xi was found by iJetha et ail d2007l) to be 
typical for moderately sized groups and poor clusters (i.e. a 
typical richness of 9 - 12). X4 an d X5 were measured by 
I Jones & Forman] d!984l) . who found them to be typical for 
moderately sized and massive clusters, respectively. Thes e 
last two models were used by IKaiser & Alexander] dT999b). 
These five estimates span three orders of magnitude and the 
full range of plausible density parameters for the host clus- 
ters of the FR lis in our sample. We expect, however, that the 
IJetha et all (120071) density estimate most closely approximates 
the ICM density of most of the host clusters in our sample 
since the clusters in our sample have a med ian richness which 
is comparable to the median richness of the lJetha et al.l(l2007l) 
sample. 3 We therefore take the \3 density model to be the 
default density model throughout this paper. 

3.3. The Jet Power Distribution 

In addition to an ICM model, the KDA model requires 
a jet power to calculate the length and luminosity as a 
function of age. To model the observed length distribu- 
tion , we therefore req uire a plausible distribution of jet pow- 
ers. IBird et afl (120081) use d two power-law distrib utions, one 
from Blun deTTet al.l ([1999) and one derived from Sadl er et aD 
(2002). Although the two distributions have very different 
slopes (0.62 and 2.6), they have similar median powers. Be- 
cause the FR II length distribution depends primarily on the 
median po wer of the FR lis , not powers at the tails of the 
distribution. iBird et al.l (120081) found similar results with both 
distributions. We find both of these distributions to be unsuit- 
able, however. If we use the KDA model to model a jet at 



3 Although Jetha et al. 1 2007) do not cite richness estimates, compiling 
richnesses for clusters in their sample from various sources in the literature 
yields a median richness of nine. 
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TABLE 3 
ICM Density Profiles 



Model a p par Environment Reference 

(kpc) (kg m -3 ) (kpc L9 kg nT 3 ) 





10 


1.7 x lCT 


-23 


1.3 x 10" 


-21 


Field & poor groups 


Blundell et al. (19991: Garrington <S 


: Conwav (1991) 


Xi 


2 


7.2 x lCT 


-22 


2.7 X 10" 


-21 


Field & poor groups 


Kaiser et al. (1997); Canizares 


et al. ( 1987) 


X3 


391 


7.2 x 10" 


-26 


6.1 x 10" 


-21 


Large groups & poor clusters 


Jetha et al. (2007) 




X4 


30 


5.0 x 10" 


-23 


3.2 x 10" 


-20 


Moderate clusters 


Kaiser & Alexander (1999b); Jones 


& Forman (1984) 


X5 


260 


5.9 x 10" 


-24 


2.3 x 10" 


-19 


Large clusters 


Kaiser & Alexander (1999b); Jones . 


& Forman (1984) 



NOTE. — Density profiles given in terms of increasing density parameter, ap@ . The density and scale length are not independent parameters 
in the KDA model, but only appeal' in their combined form ap@ . 



the median power of these distributions (^10 38 W) and with a 
length typical for the jets in our sample (^100 kpc), the pre- 
dicted radio luminosity is two orders of magnitude larger than 
the t ypical radio lumin osity of the FR lis in our sample. 

As lBird et al.l d2008) showed, lifetime estimates scale as the 
cube root of the median jet power. Nevertheless, an error in 
the power distribution by over an order of magnitude can pro- 
duce errors in the lifetime calculations by a factor of at least 
three. To accurately determine the jet lifetime, we therefore 
need a reasonable estimate of the jet powers in our sample. 

Since the KDA model is a one-to-one function of age 
and power to length and radio luminosity over the range of 
physically-motivated choices of the input and auxiliary pa- 
rameters, it is possible to numerically invert it to produce an 
age and power given a length and luminosity. For any given 
jet in our sample, we compute the set of powers which repro- 
duces the observed length over a range of ages from 10 6 yr 
to about 5 x 10 9 yr (at which point numerical calculations be- 
gin to stop converging 4 ). This upper bound is well above nu- 
merous other con straints on the lifetime of AGN, and FR lis 
specifically ("e.g.. lMartiirill2004t lO'Dea et al.ll2009l) . Because 
the observed projected length is only a lower bound on the 
true physical length of the FR II, we make an average correc- 
tion of a factor of 4/V to the length to account for this pro- 
jection effect. We similarly compute the set of powers which 
reproduces the observed luminosity over the same range of 
ages. The intersection of the loci of the two sets is then an 
estimate of the age and power of the FR II. Figure [5] illus- 
trates the locus of constant length and the locus of constant 
luminosity for a typical FR II in our sample. 

We perform this process for every jet in our sample with 
each of the five ICM models. The median powers and ages 
are presented in Table |4] and the distribution o f the sam- 
ple in power-age space for the iJetha et all (|2007) ICM den- 
sity estimate is displayed in Figure [6] These median powers 
range from 4 x 1 36 W to 10 37 W, lower than the median jet 
power assumed in B ird et all (|2008) by over an order of mag- 
nitude. These jet powers are also lower than t he minimum 
FR II jet power derived in iKaiser & Bestl d2007l) : we provide 
a plausible explanation for why their calculation is not in- 
consistent with our results in $3.21 Our results are consistent 
with other, s impler methods of estimating the jet power (e.g., 
lPunslvll2005l) . 

The resulting distribution of FR II jet powers in our sample 
is not well described by a truncated power law. To account 
for selection effects, we generate a set of mock catalogs (de- 



TABLE4 

Estimated Median Jet Powers and 
Ages 



JA^fVl UcliSlLy 


A 
r\ 


n 


t 


Model 




/ i n35 
(10 W ) 


(10 7 yr) 


(1) 


(2) 


(3) 


(4) 


V i 
A 1 


2 


211.3 


9.8 


A. 1 


3 


165.6 


7.7 


V i 
A 1 


4 


145.0 


6.6 


XI 


6 


124.7 


5.5 


Xl 


8 


113.9 


4.8 


XI 


12 


102.0 


4.0 


XI 


16 


95.1 


3.5 


X2 


2 


174.6 


13.2 


X2 


3 


134.6 


10.4 


X2 


4 


116.5 


9.1 


X2 


6 


99.0 


7.3 


X2 


8 


89.8 


6.6 


X2 


12 


79.7 


5.4 


X2 


16 


74.0 


4.8 


X3 


2 


142.9 


18.5 


X3 


3 


108.8 


14.9 


Xi 


4 


93.3 


12.8 


X3 


6 


78.5 


10.6 


X3 


8 


70.7 


9.3 


X3 


12 


62.3 


7.7 


X3 


16 


57.5 


6.9 


X4 


2 


99.4 


36.3 


X4 


3 


75.0 


29.1 


X4 


4 


63.9 


25.4 


X4 


6 


53.1 


20.8 


X4 


8 


47.6 


18.3 


X4 


12 


41.6 


15.2 


X4 


16 


38.3 


13.7 


X5 


2 


68.4 


79.2 


X5 


3 


51.9 


62.7 


X5 


4 


44.5 


54.6 


X5 


6 


37.3 


45.7 


X5 


8 


33.6 


39.7 


X5 


12 


29.7 


32.9 


X5 


16 


27.5 


29.2 



NOTE. — The median powers and ages of 
the FR lis in the full sample esti mated us- 
ing the method described in 33-31 Our de- 
fault choice of density model and axial ratio 
is shown in bold. (1): The density model (de- 
finedin Table f3). (2): The axial ratio (defined 
in 34.2) . (3): The estimated jet power. (4): 
The estimated current age. Due to selection 
effects there is no simple way to convert the 
median current age to a lifetime estimate. 



4 Specifically, for sufficiently extreme parameters, Equation (A26) of 
"HI BOOl) 1 



IKaiser & Besl I 



) fails to have a solution. 
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10 8 
Age (yr) 

FIG. 5. — A graphical explanation of the procedure to find the age and 
power of an FR II source given its length and luminosity. The solid line 
is the locus of points in age-power space consistent with the jet's observed 
luminosity, in this case 1.3 X 10 25 W Hz" 1 . The dashed line is the locus of 
points consistent with the jet's observed length, in this case 98.1 kpc. The 
intersection of the two lines is the only age and power consistent with the 
observed length and luminosity of the jet within the KDA model. The data 
shown in this plot are for a typical jet in our sample for an axial ratio of 
A = 4 and an ICM with p = 7.2 X 1(T 26 kg nT 3 and a = 391 kpc. Because 
the observed length of an FR II is shorter than its true physical length due to 
projection effects, we make an average correction to the observed jet length 
before estimating its age and power. 



scribed in further detail in §13.51 consisting of 10 6 FR lis given 
random ages, redshifts, and orientations. If we assume that 
the jet powers are distributed as a truncated power law with a 
minimum power Q m \ n and a slope a, a maximum likelihood 
analysis (described in further detail in § 14. 11 where it is applied 
to the length distribution and lifetime measurement) applied 
to the resulting power distributions i ndicates tha t the b est fit 
is gmin = 7.6 x 10 35 W, a = 2.2 for the lJetha eTail d2007b ICM 
model and an axial ratio of 4. Although better fits exist in 
principle, the fraction of detectable FR lis in them is so low 
that our sample size would be much less than 151 even if ev- 
ery BCG in the MaxBCG catalog hosted an FR II. Even in 
this best fit, a sample size of 151 FR lis can only be recovered 
if -50% of the BCGs in the MaxBCG catalog host FR lis. 
We find that the power distribution is much better described 
by a log-normal distribution. A similar analysis shows that 
the best fit has a median power of 1.1 x 10 37 W and a coeffi- 
cient of variation of 2.2 under the same density and axial ratio 
assumptions. (We perform this fit for every choice of ICM 
density and axial ratio we examine in this work. See § |4.2| and 
§14.31 for these choices.) A log-normal distribution not only 
fits the observed jet power distribution better than a power 
law, but it also does not require a very low detection probabil- 
ity, and hence does not require that a large fraction of BCGs 
harbor FR lis (—2% for a log-normal distribution vs. —60% 
for a power law). The cumulative distribution of jet powers 
is presented in Figure [7] along with the best-fitting power law 
and log-normal distributions. 

The true distribution is likely more complicated (and pos- 
sibly an approximate power law) than log-normal, but only 
appears log-normal as a result of selection effects. Although 
the mock catalogs account for any biases due to our selection 
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FIG. 6. — The distribution of the full sample (filled and open circles) and 
the volume-limited sample (filled circles only) in estimated jet power and age 
assuming an axial ratio of 4 and the X3 density model suitable for typical 
clusters (black points), the \ l density model suitable for galaxy groups (left 
group of light gray points), and the \5 density model suitable for the rich- 
est clusters (right group of light gray points). (See Table [3] for definitions 
of the density models.) Because the age and power estimates indicated here 
are determined simultaneously, errors in the power estimate can induce er- 
rors in the age estimate, and vice versa. Furthermore, although an average 
correction is made for projection effects, the age and power estimates in this 
figure do not reflect any distribution in projection angle. The age estimates 
in this figure are therefore only approximate. We use the length distribution 
to constrain the age more precisely because it is more strongly dependent on 
age than luminosity. Because the uncertainty in our estimates is dominated 
by systematic uncertainties in the ICM density and axial ratios of the FR lis, 
we do not show error bars in this plot. Systematic uncertainties are on the 
order of a factor of two in both age and jet power based on variation in these 
estimates over a reasonable range of ICM densities and axial ratios. 

algorithm, they do not account for selection effects imposed 
by the MaxBCG catalog. Specifically, the distribution of BCG 
stellar luminosities in the MaxBCG catalog can be much bet- 
ter described by a log-normal distribution than a power law. 
Since BCG stellar lumin osities are correlated with FR II jet 
powers (as we show in §15. 1.1b , it is not surprising that our 
observed distribution of jet powers resembles a log-normal 
distribution. The log-normal distribution in BCG luminosi- 
ties may be the result of the steep galaxy luminosity function 
at high luminosities and the minimum richness cutoff of the 
sample. A study of jet powers in lower-luminosity galaxies 
and field galaxies could determine if the intrinsic jet power 
distribution is log-normal for all galaxies. 

While the method described in this section produces esti- 
mates of the ages of the FR lis, we do not expect them to 
be particularly precise. This is because the estimates of the 
lengths and ages of the FR lis place equal weight on the equa- 
tion for luminosity evolution and the equation for length evo- 
lution. The length evolution of the FR II is relatively straight- 
forward to model, because it only relies on dimensional anal- 
ysis, and hence is consistent across many models. The lumi- 
nosity evolution, on the other hand, is much more difficult to 
model and so has greater uncertainty. Even if the equation 
for luminosity evolution were known perfectly, there is also a 
greater uncertainty in the flux measurement than in the length 
measurement. This is because the sources are extended and 
therefore flux from diffuse components can be lost and flux 
from spatially coincident sources can be added. We therefore 
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FIG. 7. — The cumulative distribution of jet power for FR lis in our sample 
(solid line), the best-fitting power law distribution (dotted line), and the best- 
fitting log-normal distribution (dashed line). These distributions are the ob- 
served distributions from mock catalogs which take selection effects into ac- 
count to make them directly comparable to the observed sample. The power 
law distribution shown here makes the implausible requirement that ~60% 
of BCGs in the MaxBCG harbor FR lis, but only ~7% are detected. Al- 
though better-fitting power law distributions exist, they require even smaller 
detection fractions; even if all BCGs in the MaxBCG catalog harbored FR lis, 
these distributions predict FR II sample sizes much less than the 151 we ob- 
serve. The log-normal distribution, by contrast, predicts a detection fraction 
of -60%. 

only use this process to estimate the jet power to within a fac- 
tor of a few and then use more precise methods to measure 
the lifetime. Because the lifetime calculation depends only 
on the cube root of the jet power, the uncertainty in the power 
estimate will not introduce an error in the lifetime calculation 
of more than a factor of two. 

3.4. The Q — P\ .4 Relation 

Several authors have sought a simple power- law relation- 
ship b e tween the jet p o wer and radio lumino s ity (IWillott et alj 



1999; Punsly 2005t iBirzan et all l200l ICavagnolo et"all 
20101) . Although the jet power is the fundamental physical 



parameter of interest in radio galaxies, it is not easily mea- 
surable. A simple, approximate conversion between jet power 
and the observed radio luminosity would be useful to estimate 
jet powers in large samples of radio galaxies. 

We fit the powers estimated in jj3.3| to a pow er-law in 1.4 
GHz radio luminosity. For the[jetha et all (120071) ICM density 
model and an axial ratio of 4, we find the best fit to be 



logg = 0.95(±0.03)logP L4 + 13.4(±1.1) 



(7) 



where both Q and Pi 4 are in Watts and the quoted errors are 
purely statistical. Variation of the ICM density and axial ra- 
tio changes the normalization, but leaves the exponent of the 
power law unchanged to within the quoted error. This fit is 
shown in Figure[8] The scatter in this relation is largely due to 
the fact that variations in length (and hence age) are ignored. 

This relation is st eeper than the relatio ns found by 
IBirzan et all d2008l) and lCavagnolo et al l (1201 Ol) , who find ex- 
ponents of 0.35 ± 0.07 and 0.75 ± 0.14, respectively. Our 
relation is, however, consistent with the relation between Q 
and th e rad io luminosity at 151 MHz derived by Will ott et al.l 
d 19991) and lPunslvl d2005h . who found an exponent of 0.86. 
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FIG. 8. — The best-fit power law relation between radio luminosity at 1.4 
GHz and the jet power assuming the \3 density model (defined in Table [3} 
and an axial ratio of 4. Changes to the ICM density and axial ratio change 
the overall normalization but do not have a large effect on th e exponent. The 
method for calculating the jet powers is described in j)3.3l We emphasize 
that the best fit relation is heavily model-dependent due to our use of the 
KDA model to convert from jet length and radio luminosity to the jet power. 
This relation is therefore best interpreted more as a predic tion of the KDA 
model than a direct measurement of the Q — Pi a relation. 

We emphasize that our fit is heavily model-dependent be- 
cause our conversion from jet length and radio luminosity to 
jet power is derived from the KDA model. While the KDA 
model likely provides estimates of jet powers which are cor- 
rect to within a factor of a few, any errors in the relationship 
between jet power and radio luminosity derived in the model 
will be reflected in the exponent of our fit. The scatter in this 
relationship will also be influenced by the intrinsic jet power 
and age distribution of the sample. The scatter for our sample 
may be relatively modest because the jet power distribution 
follows a log-normal distribution, which may not be the case 
for all FR II samples due to the selection of the input cluster 
catalog (see j32.lt . 

3.5. Mock Catalog Generation 

We compare the observed length distribution of FR lis to 
mock catalogs to infer their intrinsic properties. These mock 
catalogs are based on large, simulated datasets that sample a 
suite of KDA model parameters (e.g., the ICM density and 
axial ratio) and a range of maximum ages. For any given 
choice of KDA model parameters, we create 3 1 mock cata- 
logs with maximum ages ranging from 1.1 x 10 7 to 1.2 x 10 9 
yr and distributed evenly logarithmically. Roughly 1.3 x 10 6 
jets are then created by randomly choosing an age, power, and 
redshift for each jet. The age is sampled from a uniform dis- 
tribution ranging from to the maximum age chosen for that 
particular catalog and the jet power is sampled from the log- 
normal distribution derived in 93.31 To account for selection 
effects, we decrease the median power by 15%, an empiri- 
cally determined amount for which the jet power distribution 
of the mock catalogs best matches the observed jet power dis- 
tribution. The redshift of the simulated jet is then sampled 
from the redshift distribution of the entire MaxBCG catalog. 
(The redshift distribution is necessary when calculating the jet 
luminosities due to the effect of the CMB energy density on 
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synchrotron losses.) We assume that the birth rate of FR lis 
is constant across the redshift range of the Max BCG catalog. 
The support for this assumption is discussed in 35.31 

Once a jet has been created, we calculate its .^-corrected lu- 
minosity and randomly orient it on the sky to calculate its pro- 
jec ted l ength. We then apply the selection criteria described 
in 32.3l to create a mock catalog for that parameter set. 

4. THE FR II LIFETIME 

4.1. Maximum Likelihood Fitting 

Our approach to estimate the maximum lifetime for FR lis 
is to compare the projected length distribution of the mock 
catalogs to the observed distribution. We fit the observed dis- 
tribution to the mock catalogs with the maximum likelihood 
method. The result is the best-fit maximum age for a given set 
of model parameters. 

The likelihood of any maximum age is determined by 
constructing a probability distribution function (PDF) of the 
length distribution for a mock catalog with 1-kpc-wide bins. 
If the number of jets in the mock catalog that fall into the i* 
bin is n,, the probability that a jet will fall in that bin is rii/N, 
where N is the total number of jets in the mock catalog, £\ n,-. 
With a PDF, /, so defined, the likelihood of any model is then 

C = l[f(lj), (8) 

j 

where lj is the projected length of the / h jet in the observed 
sample. For any given choice of parameters, we take the best- 
fit lifetime to be the maximum age of the mock catalog which 
maximizes the likelihood. The error on the lifetime calcula- 
tion is taken to be the FWHM of the likelihood distribution. 

Figure [9] shows likelihood distributions for several mod- 
els. The length distribution of the mock catalog for a sam- 
ple model at three ages is displayed with the observed length 
distribution in Figure [10] Because the projected length dis- 
tribution is extremely sensitive to age, we can determine the 
FR II lifetime with relatively high precision. Table [5]presents 
the lifetimes and their uncertainties for every model used. We 
note that we only measure the typical FR II lifetime and there 
is likely some intrinsic dispersion which is larger than the sta- 
tistical uncertainties indicate. We discuss any correlations be- 
tween the typical FR II lifetime and various jet properties in 
35.1.41 The best-fitting lifetime for our default ICM model is 
1.9xl0 8 yr. 

4.2. Dependence on Axial Ratio 

We define the axial ratio as the ratio of the length of the 
lobe to the perpendicular distance between the edge of the 
lobe and its axis. (Some other studies define the axial ratio 
as the ratio between the length of the lobe and the width of 
the lobe; our definition is twice this other convention.) Due 
to projection effects, background noise, and the fact that the 
shapes of real radio lobes are more complicated than the sim- 
ple geometry of radio lobe models, unbiased and robust mea- 
surements of the axial ratios of FR II sources are difficult to 
obt ain. Nevertheless, nearly all estimates fall betwe en 2 and 
16 dLeahv & Williamsl [l984t iMachalski et alJl2004l) . Some 
studies of extremely powerful jets have found even larger ax- 
ial ratios, but there appears to be a relationship between jet 
power and axial ratio, with more powerf ul jets having unusu- 
ally large axial ratios dLeahv et al. 10*989 ). While it is possible 
to estimate the axial ratios of the FR lis in our sample from the 
FIRST images, the width of the lobes are in many cases com- 
parable to the size of the FIRST PSF, so such estimates only 
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provide lower bounds on the axial ratios. Since our sample 
does not include atypically powerful jets, we simply examine 
axial ratios within the entire plausible range of 2 < A < 16 and 
note that our estimates are most consistent with an axial ratio 
of 3 - 4, but could be larger. Higher-resolution data of some 
of the FR lis in our sample would be needed to accurately 
measure their axial ratios. 

To first order, one expects that as the axial ratio increases, 
the lifetime should decrease. At high axial ratios, the jet does 
not need to expend as much energy inflating a fat cocoon 
around it and can therefore drill through the ICM more ef- 
ficiently; thus, a jet of given power with a large axial ratio can 
reach the same length as a jet with a small axial ratio in less 
time. But this simple relationship is complicated by the fact 
that the axial ratio changes the estimated power as well. A 
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FIG. 9. — Likelihood distributions of three models for a range of FR II 
lifetimes. The peak of the likelihood distribution for a particular model is 
the lifetime that best matches the observed data. The solid points represent 
the likelihood distribution for the default ICM density model \3 an d an axial 
ratio of 4. The open circles connected by a solid line represent the default 
ICM density model X3 and an axial ratio of 12. The open circles connected 
by a dotted line represent the \4 ICM density model and an axial ratio of 4. 
See Table[5]for a description of the ICM density models. 

jet with a larger axial ratio requires less power to produce a 
given length in a given time, so larger axial ratios will lead 
to lower power estimates. This will partially counterbalance 
the increased age estimate and produce a smaller axial ratio 
dependence on age than one would naively expect from the 
jet length and age alone. Nevertheless, because the lifetime 
is relatively insensitive to the power, this correction does not 

d ominate. 

iBird et aD (120081) found that the axial ratio has a stronger ef- 
fect on FR II lifetime calculations than all other KDA model 
parameters. To evaluate the dependence of our FR II lifetime 
measurements on axial ratio, we calculated the lifetimes of 
jets in our sample for different axial ratios. For any particular 
choice of axial ratio, we used the KDA model to estimate the 
jet powers, generated a set of mock catalogs across a range of 
ages, and fit these to the observed length dis tribu tion with the 
maximum likelihood method described in jj4.ll The depen- 
dences of power and lifetime measurements on the assumed 
axial ratio are shown in FiguresQT]and[T2] The inferred power 
decreases as the axial ratio increases, but the effect is only ap- 
proximately a factor of two decrease over a factor of eight in- 
crease in axial ratio. The age estimate similarly decreases as 
the axial ratio increases, with the effect also being a roughly 
factor of two decrease over a factor of eight increase in axial 
ratio. This effect is even weaker when the change in inferred 
jet power is taken into account. 

4.3. Dependence on the ICM Model 

The model length distribution also depends on the ICM 
density parameter. ICM density measurements of group and 
cluster samples can vary by several orders of magnitude 
(Vikhl inin etai]|2006t Uetha etalj 120071: iFreeland & Wilcotsl 
1201 ll) . Since the density of the model ICM can potentially 
have a large effect on FR II lifetime estimates, we produced 
sets of mock catalogs spanning about two orders of magnitude 
in density. The density profiles we use and their references are 



listed in Table [3] 

At higher densities, a jet of a given power will require more 
time to grow to a given length, so to first order, increased den- 
sities will lead to increased lifetime estimates, approximately 
linearly. However, increased densities also increase the lumi- 
nosity of a jet at a fixed power; hence higher densities lead to 
lower power estimates. This, in turn, will result in yet longer 
lifetime estimates in second order. 

For each ICM model, the length distribution of the mock 
catalogs was fit to the observed length distribution of the 
flux-limited sample with the maximum likelihood method de- 
scribed in {34.11 The dependence of power and lifetime mea- 
surements on the assumed ICM model are shown in Fig- 
ures [13] and [14] These figures indicate that larger densities 
yield larger lifetime estimates and smaller power estimates, 
though in the case of the jet power estimates the dependence 
appears to be more complicated in detail than a power law. 

4.4. The FR II Duty Cycle 

Because our mock catalogs take into account the biases in- 
duced by selection effects on our sample, we use them to cal- 
culate a robust upper bound on the duty cycle of FR lis in 
BCGs, where the FR II duty cycle is the average FR II life- 
time divided by the fraction of time that the jet remains active. 
Although it is not possible to directly measure the fraction of 
time that a jet remains active, a lower limit on this fraction 
can be found in a straightforward manner. If we assume that 
every BCG in a cluster has a radio jet at some point in its life, 
then the fraction of time that a jet is active is the number of 
BCGs with radio jets divided by the total number of BCGs. 
Since it might be that many BCGs never have a radio jet, this 
fraction is necessarily a lower limit. 

The total number of BC Gs in our sample is 13,8 23, the size 
of the MaxBCG catalog of iKoester et alj d2007al) . The num- 
ber of BCGs with FR lis is our sample size of 15 1 FR lis plus 
some additional FR lis that fall below our selection criteria. 
To estimate the number of FR lis that are missed, we calcu- 
late the fraction of simulated jets that pass all of the detec- 
tion thresholds in our mock catalogs. The mock catalogs are 
populated with jets resembling the underlying population of 
FR lis, so the fraction of jets which pass our detection cuts is 
an estimate of the completeness of our sample. For the best-fit 
age and the default parameters, this completeness fraction is 
61.7%. The overwhelming majority of missed sources are too 
short. We divide the number of jets in our sample by this com- 
pleteness fraction and estimate that there are approximately 
245 FR lis in these BCGs. The fraction of time that a jet is 
active is then 245/13,823, or 1.8%. This frac tion is in agree- 
ment with the value found bv lBird et al.l (120081) for galaxies in 
the group environment. Combined with our best estimate of 
the FR II lifetime of 1.6 x 10 8 yr, this yields an upper bound 
for the duty cycle of 9 x 1 9 yr. As this upper bound is compa- 
rable to a Hubble time, it implies that most BCGs undergo an 
FR II phase only once. If a substantial fraction of BCGs are 
never FR lis (either being always inactive or only appearing 
as FR Is), then FR lis could be episodic. 

5. THE RELATIONSHIP BETWEEN FR IIS AND THEIR 
ENVIRONMENTS 

We here determine whether any of the properties of FR lis 
determined in the above analysis are correlated with the FR II 
environment. In ^5.11 we examine correlations of jet power 
with stellar luminosity, BCG velocity dispersion, cluster rich- 
ness, and FR II lifetime. The impact of our assumptions about 
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FIG. 10. — A comparison of the cumulative projected length distribution of the sample (solid line) to the cumulative length distributions of a model for three 
different lifetimes (dashed and dotted lines) assuming the default ICM density \3 ( see Table[5J and an axial ratio of 4. The dashed line shows the FR II lifetime 
from the model with the maximum likelihood, which in this case was 1.9 X 10 s yr. The left and right dotted lines represent FR II lifetimes of 6.1 X 10 7 yr and 
3.4 X 10 8 yr, respectively. 

the ICM model used in 35. II is assessed in 35.21 We finally 
examine the broader relationship between the FR II and radio 
galaxy fractions and galaxy cluster richness in § 15.31 

5.1. Correlations with Jet Power 

Radio jets are powered by the accretion of gas onto nu- 
clear, supermassive black holes (SMBHs). Although the 
mechanism by which infalling gas is expelled and collimated 
from the SMBH is not well understood, most jet models rely 
on some variant of the B landford-Znajak (BZ) mechanism 
(Bland ford & Znajekll 1977b . The BZ mechanism explains jet 
outflows from SMBHs by postulating a magnetic field an- 
chored to an accretion disc around a rapidly spinning SMBH. 
As the accretion disk rotates around the SMBH, the field lines 
become wound around the spin axis. Some charged particles 
falling into the SMBH will then follow the field lines and es- 
cape along the spin axis of the SMBH. Jet models which apply 
the BZ mechanism find that the power that can be extracted 



from the SMBH and injected into the jet is roughly propor- 
tional to the accretion r ate and the square of the SMBH spin 
dMcNamara et al.|[201 lb . Larger black holes should therefore 
generally produce more powerful jets. It is interesting to see 
whether we can recover this relationship in our FR II sample. 
While direct measurements of SMBH masses are not avail- 
able at the distances in our sample, we can check for the exis- 
tence of a correlation between jet power and two proxies for 
SMBH mass: BCG luminosity and BCG velocity dispersion. 
We also search for a correlation between jet power and cluster 
richness, which may be related to the fuel source and thus the 
accretion rate, and between lifetime and cluster properties. 

Because the MaxBCG catalog is not complete and our FR II 
sample is flux-limited, selection effects could induce an arti- 
ficial correlation between jet power and BCG or cluster prop- 
erties. For instance, if the MaxBCG catalog suffered from 
Malmquist bias, a distant FR II in our sample would have to 
be both luminous at radio wavelengths (and hence have a large 
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Axial ratio 

FIG . 11. — The dependence of the axial ratio on the median inferred power 
of the FR II sample. The points along each dotted line represent the median 
powers inferred for a particular density model. From bottom to top, the den- 
sity models used are xi - Xs- See Table fj] for the details of each density 
model. 
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FIG. 12. — The dependence of FR II lifetime estimates on the assumed 
axial ratio. Each dotted line represents a set of lifetime estimates for a single 
density. Lower tracks assume lower densities. The densities used are listed 
in Table [5] The solid line represents the axial ratio dependence one would 
expect if one only considered the effect of axial ratio on the length of the 
jets as a function of time. Because the axial ratio also has an effect on our 
power estimates of the jet which tends to counteract the effect of axial ratio 
on length evolution, the final lifetime estimate is less sensitive to the axial 
ratio than found by other studies (e.g., Bird et al. 2008). 

jet power) and be hosted by a luminous galaxy. This would 
lead to a concentration of FR lis with high jet power and high 
BCG optical luminosity in the sample and thereby induce a 
correlation between the two properties. Completeness tests 
based on mock catalogs indicate that the MaxBCG catalog is 
principally biased against low-mass clusters since the member 
galaxies of such clusters tend to be less significant overden- 
sities, although there is a slight bias against extremely hi gh- 
mass clusters as well (see Figure 7 of K oester et al.ll2007al for 
the detailed completeness function of the MaxBCG catalog). 
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FIG. 13. — The dependence of the FR II lifetime calculation on the density 
of the ICM. In each case the density profile is described by a power law with 
a slope of 1.9. Each track represents a different choice of axial ratio; lower 
tracks have larger ax ial r ati os. T he values of the axial ratios are the values of 
the points in Fi gures [TTl and IT2] 
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FIG. 14. — The dependence of the estimated FR II jet power on the density 
of the ICM. Each track represents a different cho ice of a xial ratio; lower 
tracks have larger axial ratios. Refer to Figures [TTl and H2l for the axial ratio 
of each track. 



Low-mass clusters typically host lower-luminosity BCGs, so 
we should expect some correlation between jet power and 
BCG luminosity due to these selection effects. To avo id th is 
bias, we use the volume-limited sample described in 32.41 to 
search for correlations between measured properties of the 
FR lis and properties of their hosts. We note that while we 
do not remove the selection biases of the MaxBCG catalog, 
selection biases must be present in both variables to induce a 
correlation; mitigating the selection biases in just one variable 
is sufficient to remove this correlation. 
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FIG. 15. — The correlation between the r-band luminosity of the host BCG 
and the esti mated jet power. Filled circles represent jets in the volume-limited 
sample (see j|2.4t and open circles represent jets in the full sample, but not the 
volume-limited sample. The correlation between the r-band luminosity and 
the jet power for jets in the volume-limited sample is 0.284; there is a 0.3% 
chance that this correlation could be drawn from an uncorrelated population. 
When the jets from the full sample are included, the correlation becomes 
stronger and more statistically significant. 

5.1.1. The Correlation between Jet Power and Stellar 
Luminosity 

The MaxBCG catalog provides a /^-corrected r-band and 
/-band luminosity for every BCG in the sample. Because el- 
liptical galaxies generally exhibit very little star formation and 
have little dust, both the r- and /-band luminosities should be 
well-correlated wi th the stellar mass , which is correlated with 
the SMBH mass jNovak et alj|2006l and references therein). 
If there is some power-law relationship between BCG lumi- 
nosity and SMBH mass, we expect a correlation between the 
logarithm of the BCG luminosity and the logarithm of the jet 
power. We cal culate that the Pearson correlation coefficient 
ilPearsonll 18951) between logg and logL r is 0.284. Given that 
our volume-limited sample has 106 degrees of freedom, the 
probability that a correlation coefficient of 0.284 or higher 
could be randomly drawn from uncorrelated data is only 0.3% 
and therefore statistically significant. A correlation of similar 
strength exists between the /-band luminosity and jet power. 
Since jet power is very strongly correlated to radio luminos- 
ity (see 33.31 >. a correlation of similar statistical significance 
exists between the radio luminosity of the radio lobes and the 
optical luminosity of the host galaxy. We show the relation- 
ship between jet power and r-band luminosity in Figure [15] 

5.1.2. The Correlation between Jet Power and BCG Velocity 

Dispersion 

The velocity dispersions of elliptical galaxies are 
tightly correlated w ith their nuclear SMBH masses (e.g., 
Gulte kin etal1l2009h . If there is a correlation between jet 
power and SMBH mass, we also expect a correlation between 
jet power and host velocity dispersion. We calculate the Pear- 
son correlation coefficient between log Q and log a and find it 
to be only 0.1 16. There is a 33% chance that this correlation 
could be drawn from an uncorrelated population, and is there- 
fore not statistically significant. This result is at odds with 



the correlation between jet power and BCG stellar luminosity 
and may be due to two reasons. First, SDSS does not have 
spectra of the BCGs of all of the FR lis in our sample. We 
only have velocity dispersions for two-thirds of the BCGs 
in the volume-limited sample, which reduces the statistical 
power of our sample. Second, and likely more important, is 
that although the stellar luminosities of elliptical galaxies are 
strongly correlated with their velocity dispersions, velocity 
dispersion is a weak function of luminosity (proportio nal to 
the fourth root of luminosity; iFaber & Jackson! Tl 976). Our 
FR II sample has a dynamic range in BCG luminosity of 
roughly an order of magnitude. Even if jet power and optical 
luminosity were perfectly correlated, the resulting correlation 
between jet power and velocity dispersion would exhibit a 
dynamic range of less than a factor of 2. This is small enough 
to be overwhelmed by moderate scatter. 

5.1.3. The Correlation between Jet Power and Cluster 
Richness 

Richer clusters tend to host larger BCGs that, in turn, tend 
to host larger SMBHs. Richer clusters also have more hot gas 
and many clusters have cooling times shorter than 1 Gyr. With 
larger SMBHs and a larger reservoir of accretion material, it 
is plausible that cluster richness would be correlated with jet 
power. Alternatively, such a correlation could be present due 
to a correlation between jet power and BCG luminosity paired 
with a correlation between BCG luminosity and cluster rich- 
ness. We calculate that the Pearson correlation coefficient 
between logg and log/V ga i is 0.094. There is a 33% proba- 
bility that this correlation could be drawn from an uncorre- 
lated population and therefore it is not statistically significant. 
The non-detection is not too surprising as we expect that this 
correlation would be weak at best since cluster richness is a 
global property of the cluster and the jet power depends on the 
microphysics governing the immediate vicinity of the SMBH. 

5.1.4. Correlations with FR II Lifetime 

The FR II lifetime is indicative of the length of time that a 
galaxy can successfully supply its SMBH with cold gas. Since 
the mechanism that supplies this gas is not well understood, 
we use our data to determine if FR II lifetime is correlated 
with any BCG or cluster properties. We test for a correla- 
tion in two different ways. The first method is identical to 
the method for measuring correlations with BCG and cluster 
properties explained abo ve. When we estimate the FR II jet 
power (described in $3.31 >. we simultaneously estimate a crude 
(due to projection effects) FR II age for every jet in our sam- 
ple. We then calculate the Pearson correlation coefficient be- 
tween these nominal ages and various BCG and cluster prop- 
erties (such as the BCG stellar luminosity). We find that none 
of the correlations between age and the BCG and cluster prop- 
erties discussed above are statistically significant. The major 
drawback of this method is that the age estimate is sensitive 
to the length of the FR II and, due to projection effects, the 
measured projected length can be substantially smaller than 
the lobe length. While we make average corrections for pro- 
jection effects, they nevertheless introduce scatter into the re- 
lationship between the FR II age and the various BCG and 
cluster properties. Furthermore, any individual FR II is al- 
most equally likely to be detected at any point in its lifetime. 
Since we use the FR IPs estimated age (not its lifetime), this 
will also introduce significant scatter. 

Our second approach is to split the sample into two bins at 
the median of each parameter and test for correlations. This 
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provides two length distributions for each parameter. We per- 
form the maximum likelihood analysis described in $4. H and 
fit these distributions to the length distributions of the mock 
catalogs and see if there is a statistically significant difference 
in lifetimes between FR lis in the two bins. However, we 
find no difference in lifetime between the two bins when we 
split by stellar luminosity, BCG velocity dispersion, or cluster 
richness. 

5.2. The Impact of the ICM 

The above analyses were all performed under the straight- 
forward but unrealistic assumption that the ICM is identical 
for all clusters in our sample. It is known, however, that richer 
clusters tend to harbor denser ICMs with larger scale radii. 
We therefore split our sample into two bins at the median rich- 
ness (14) to reexamine how strongly the FR II jet power and 
lifetime relate to the host cluster ICM. The low richness bin 
has richness values from 10 to 14. We assign th ese clusters 
the int ermediate ICM model (x3 m Table [5) from lJetha et al.l 
(2007), which is derived from large groups and small clusters 
with a median richness of nine. Richer clusters (greater than 
14) are assigned the densest model (xs in Table [3), which is 
the geo metric mean of the densities found by Jones & Forman 
( 1984) for clusters with a median richness of 20. 

For each FR II we estimate the jet power with the method 
described in 33. 31 and the density appropriate to the richness 
of the cluster. The power distribution as a function of richness 
is shown in Figure[l6] Although there is a correlation between 
richness and FR II jet power that is statistically significant at 
the two-sigma level, it is weak. If we simply compare the 
mean power of all jets less than or equal to the median rich- 
ness with those greater than the median richness, we find that 
the jets in richer clusters are, on average, less powerful by a 
factor of 2. Although this density parameterization is some- 
what artificial because the ICM density is expected to steadily 
continue to increase with richness, albeit with substantial scat- 
ter, and spans an extreme range in density, this parameteriza- 
tion should be sufficient to identify a strong correlation. As 
the correlation between inferred jet power and cluster richness 
is only marginal, both this result and those in 35.1 l imply that 
the properties of FR lis are relatively insensitive to large-scale 
environmental factors. 

5.3. FR II and Radio Galaxy Fraction vs. Cluster Richness 

If radio-mode feedback can prevent substantial gas cooling 
in clusters, there may be an increase in duty cycle, lifetime, 
and jet power with cluster richness. In the previous subsec- 
tions, we showed there is no strong correlation between jet 
power or lifetime and richness. Here we investigate differ- 
ences in duty cycle with the fraction of clusters with an FR II 
source and the fraction of clusters with a radio source of any 
kind. We calculate the FR II fraction in richness bins that have 
been corrected for undetected sources. This process would 
eliminate any potential spurious correlations between FR II 
fraction and cluster richness due to the correlation between 
cluster richness and redshift. 

To calculate the fraction of all radio sources we cross- 
correlated the FIRST and MaxBCG catalogs. Rather than 
select radio sources at least 10 "8 (two FIRST resolution el- 
ements) from the BCG as for FR lis, we selected BCGs with 
a radio source within 10 "8. To mitigate Malmquist bias, we 
furthermore selected only those radio sources whose radio 
luminosity was greater than 1.7 x 10 22 W Hz" 1 at 1.4 GHz. 




Richness 

FIG. 16. — The FR II mean jet power as a function of host cluster richness 
after applying a parameterization between cluster richness and ICM density. 
Since larger clusters tend to have denser ICMs, we split our sample along the 
median cluster richness of 14 and assume the X3 ICM model (see Table [5) 
from Jetha et al. (2007) for the low-richness clusters (open circles) and the 
Xs ICM model from lones & Forman ( 1984) for the high-richness clusters 
(filled circles). We then estim ate the jet power of each FR II with the inferred 
density as described in j)3.3l The resulting correlation between jet power and 
cluster richness, while statistically significant, is weak. 

This corresponds to the lowest luminosity that FIRST can de- 
tect at z = 0.3, th e largest redshift in the MaxBCG catalog. 
iBird et al. (2008) performed a similar analysis using two ad- 
ditional, more stringent luminosity cutoffs and found no sub- 
stantive difference in the relationship between radio fraction 
and cluster richness or luminosity. We therefore only use this 
luminosity cutoff. 

The FR II and radio fractions as a function of cluster rich- 
ness are shown in Figure [T7] There is an increase in radio 
fraction by a factor of 1.5 - 2 with richness. Although the 
statistics on the FR II fraction are weaker due to its smaller 
sample size, there appears to be a corresponding trend in FR II 
fraction with cluster richness. We perform a similar analy- 
sis to examine the relationship between the FR II and radio 
fractions with BCG luminosity and find that both fractions in- 
crease strongly with BCG luminosity. These relationships are 
shown in Figure [18] The increase in the radio fraction with 
stellar luminosity is similar to t he increase in the radio frac- 
tion with stellar mass shown bv lBest et al.l (|2007) and ?. We 
also find that the FR II fraction does not evolve over the range 
of redshifts in the MaxBCG catalog. We show the FR II frac- 
tion across the redshift range of the volume-limited sample in 
Figure [19] 

6. DISCUSSION 

6.1. Comparison to FR II Lifetime Measurements in the 
Literature 

Our measurement of the FR II lifetime is mo st directly 
compa rable to the FR II lifetime measurement by IBird et all 
(2008) as their approach is largely identical to our own. They 
also modeled the observed FR II length distribution with 
mock catalogs to account for selection and projection effects 
and used the KD A model to determ ine the length and lumi- 
nosity evolution. IBird et all ([2008) found an FR II lifetime 
of 1.5 xlO 7 yr, which is an order of magnitude below our 
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FIG. 17. — The FR II fraction (filled circles) and fraction of all radio 
sources (open circles) as a function of cluster richness in the MaxBCG cat- 
alog. We define radio sources to be BCGs with radio emission detected by 
FIRST within 1()"8 (i.e. two FIRST resolution elements). The horizontal er- 
ror bars are the bin widths and the vertical error bars are the binomial errors. 



0.04 



0.03 



^0.02 



0.01 




FIG. 18. — The FR II fraction (filled circles) and radio fraction (open cir- 
cles) as a fu ncti on of the BCG r-band luminosity. Radio sources are defined 
as in Figure [IT] The horizontal error bars are the bin widths and the vertical 
error bars are the binomial errors. 

value of 1.9 xlO 8 yr. T here are three major differences be- 
tween iBird et all (120081) and this work that togeth er account 
for thi s difference. The first is that the FR lis in the lBird et al.1 
(2008) sample were selected from a catalog of galaxy groups 
rather than gala xy clusters. Since the host group richness was 
much smaller, IBird et al.1 (120081) took the \i ICM model as 
their default rather than the X3 model, which has a factor 
of two lar ger density param eter. The second major differ- 
ence is that Bir d et all d2008l) assumed jet power distributions 
drawn from the literature rather than deri ving the jet power 
distribution from their (smaller) sample. IBird et al.l ( 2008 ) 
used jet power dis tributions from Blundell et al. (19991) and 
Sadler et al] (120021) and found that the lifetime is insensitive 
to these two choices. This is because the lifetime calcula- 
tion depends more strongly on the median value of the jet 
power distribution than on the higher-order moments; since 
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FIG. 19. — The FR II fraction of the volume-limited sample across the red- 
shift range of the MaxBCG catalog. Horizontal error bars reflect bin widths 
and vertical error bars reflect the binomial error on the FR II fraction in each 
bin. The mean FR II fraction across the entire sample is shown by the dot- 
ted line with its binomial error represented by the shaded region. There is 
no evolution in the FR II fraction over the redshift range of the MaxBCG 
catalog. 

the Blu ndell and Sadler distributions have similar median 
powers. IBird et al.l d2008) point out that it is unsurprising that 
the two distributions produce similar lifetimes despite their 
different shapes. However, these distributions were derived 
from sam ples of mu c h mor e luminous FR lis than those in 
either the IBird et al.l (120081) sample or our own; FR lis in 
iBlundell et al J (119991) were selected from the 3C, 6C, and 7C 
catalogs and the ISadler et al.l d2002l) distribution was derived 
from the bright e nd of the luminosity function. Consequently, 
IBird et aD (120081) overestimated the median jet powers of the 
FR lis in their sample by approximately two orders of mag- 
nitude, and this led to an underestimate of the lifetime by a 
factor of approximately five. The final difference between the 
two studies is our use of an updated coef ficien t c\ to calcu- 
late the lobe length (described in detail in $3. 11 1. The updated 
coefficient is sm aller than the original coefficient c\ used by 
IBird et all (120081) by about 40% and this increases the life- 
time estimates by about 70%. When the differences between 
the two studies in jet power distribution and c\ are taken into 
account, they almost entirely explain the order-of-magnitude 
discrepancy in the lifetime measurement. The remaining dif- 
ference between the two measurements is less than a factor of 
two and can be due to uncertainties in the two measurements, 
or a small intrinsic difference in the FR II lifetime from group 
to cluster environments , or both. 

IBlundell et al.l d 19991) estimated the maximum lifetime of 
FR lis from their observed distribution in the P-D plane. They 
noted that if the lifetime were too short, one would predict far 
fewer long FR lis than are observed and there would be a 
sharp maximum-size cutoff. However, if the maximum life- 
time were too large, there would be no effect on the distribu- 
tion of points in the P-D plane for a flux-limited sample since 
FR lis above the true maximum lifetime would generally fall 
below the flux limit; larger maxi mum lifetimes only de crease 
the predicted detection fraction. IBlundell et al.1 ([1999) found 
that the smallest maximum lifetime consistent with their ob- 
servations was 5 x 10 8 yr, in reasonable agreement with our 
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value. The similarities in these lifetime estimates are interest- 
ing because our study probes FR lis with jet powers approxi- 
mately two orders of magnitude below those in the sample of 
iBlundell et aT] d!999h . IWang & Kaisel d2008l) fit a very simi- 
lar sample with a typical maximum lifetime of a few x 10 7 yr , 
although that choice was motivated by the iBird et al.1 (2008) 
v alue. 

lO'Dea et al.1 d2009b found FR II lifetimes of a few x 
10 6 ~ 7 yr with the spectral index gradient method dLeahy et al.l 
1989). They used a sample of 31 FR lis at the high end of the 
FR II luminosity function (all of the FR lis in their sample are 
in the 3C catalog), whereas sources in our sample are typically 
at much lower luminosities. As they note, because their flux- 
limited sample only includes the most powerful sources at 
high redshift, these sources are expected to have the smallest 
lifeti mes. The observed e xpansion velocities of the sources 
in the lO'Dea et al.l J2009) sample are larger than the expan- 
sion velocities observed in our sample by about two orders of 
magnitude. This difference does not appear to be due to differ- 
ences in the two methods, however, since we obtain expansion 
velocities consistent with those found bv lO'Dea et alj d2009l) 
when we apply the KDA model to the sources in their sample. 
This suggests that the discre pancy in the lifeti me measure- 
ments between our work and lO'Dea et al.l(l2009l) is due to the 
difference in jet powers between the two samples rather than 
systematic discr epancies between t he two methods. 

Additionally, IWang et all (1201 ll) found that the length of 
FR lis is limited to a maximum value b y the entrainment of 
gas within the lobe. I Wang et all (1201 ll) argued that beyond 
this maximum length the radio jet is disrupted, transforming 
the FR II into an FR I. Because more powerful jets reach this 
maximum length more quickly, lifetime studies at the high 
end of the FR II luminosity function may systematically find 
lower lifetimes. 

6.2. FR lis as Sources of Heating in the ICM 

An outstanding problem in high-mass galaxy formation is 
the suppression of substantial gas cooling. As BCGs are al- 
most always at the centers of galaxy clusters, cool gas from 
the surrounding ICM is predicted to condense onto the BCG. 
Were this process to proceed uninterrupted, BCGs would 
grow to be much larger than observed. To quench these cool- 
ing flows, some mechanism must exist to heat the ICM. 

To first approximation, the heating of the ICM due to a radio 
lobe is simply the enthalpy of the system, 
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where 7 is the adiabatic index, V is the volume of the lobe, p 
is the pressure of the ICM imme diately outside the l obe, and f 
is an ad hoc uncertainty factor dBirzan et al.ll2004t iBest et all 
2007). Although we do not know the pressure of the ICM 
outside the lobe, the KDA model does provide the pressure 
inside the lobe (equation A4 of Kaiser & Best, 2007). Since it 
is reasonable to assume that the lobe is overpressured relative 
to the ICM, we calculate the enthalpy of the system with the 
lobe pressure and thereby obtain an upper limit. 

To calculate the volume of the lobe, we assume a cylindri- 
cal geometry with length / and an axial ratio of A = 4. The 
observed lengths of the FR lis are all extended by a factor of 
7r/4 to correct for projection effects. Finally, since the gas is 
relativistic, 7 = 4/3. The enthalpy is therefore 



We use the KDA model to calculate the lobe pressure for 
every FR II in our sample based on the estimated ages and 
powers of the jets (see 33.3b . Although we may derive a lower 
limit on the total ICM heating due to the adiabatic expansion 
of the lobe by dividing the enthalpy by the estimated age of the 
FR II, such an estimate is only partially useful in determining 
whether FR lis can quench cooling flows in galaxy clusters. 
To be a successful mechanism to quench cooling flows, FR 
lis must not only provide a sufficient amount of total heat, but 
must also provide a sufficient amount of heat throughout the 
central region of the cluster. Otherwise, cooling flows would 
be quenched in some regions and would proceed unabated in 
others. In particular, if gas on the equatorial plane between 
the two lobes is not shocked, the assumption that the lobe 
expands adiabaticly will provide an upper limit on the heating 
of gas in this region. 

To determine whether ICM heating from the adiabatic ex- 
pansion of the lobe is sufficient to halt substantial gas cooling, 
we compare the average heating rate to the average cooling of 
the ICM. While we do not have X-ray luminosities to directly 
estimate cooling rates, we can estimate the typical cooling 
rate from the scaling relation between X-ray luminosity and 
cluster velocity dispersion of 



log 10 (L 45 ) = -1.34 + 2.0 x log 10 ((T5oo) erg s 
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where L45 is the X-ray luminosity in units of 10 45 erg s" 1 and 
0500 is the cluster ve locity dispersion in units of 500 km s" 1 
dOrtiz-Gil et al.ll2004b . We estimate the cluster velocity dis- 
persio ns from the scaling relation measured by iBecker et al.l 
d2007l) of 

(lncr) =6.17 + 0.4361n/V g 2 S°/25. (12) 

These two equations provide an estimate of the ICM cool- 
ing in the MaxBCG clusters. The heating rate from the en- 
thalpy of the lobes (assuming / = 1) is shown with the esti- 
mated cooling rate in Figure [20] The typical heating rate due 
to the FR II lobes is at least an order of magnitude smaller 
than the estimated cooling rate, implying that heating from 
radio-mode feedback from FR lis is insufficient to quench 
cooling flows in rich clusters. Figure [20] also demonstrates 
that the total jet powers exceed the amount required to coun- 
teract cooling by approximately an order of magnitude. That 
is, only ^10% of the jet power is needed to counteract the 
predicted cooling, while only ~1% of the jet power appears 
to contribute to the enthalpy of the lobes. These jets are con- 
sequently sufficiently powerful to counteract cooling in indi- 
vidual clusters, but this power may not couple to the ICM. 

There are several sources of uncertainty in the enthalpy cal- 
culation that may bring the heating and cooling rates into bet- 
ter agreement. First, our calculation of the enthalpy of the 
lobe requires knowledge of the axial ratio of the lobe. We 
assume an axial ratio of A = 4, but because H oc A" 2 , an over- 
estimate of the axial ratio will lead to an underestimate of the 
enthalpy. However, the axial ratio can generally be no smaller 
than A = 2, which corresponds to spherical lobes. This would 
increase the enthalpy by a factor of four. 

Second, our calculation of the enthalpy of the lobe requires 
knowledge of the pressure within the lobe. This can be calcu- 
lated from the KDA model, but is dependent upon the density 
of the ICM outside the lobe. Specifically, if the density profile 
of the ICM has a power-law slope of (3 = 2, p oc pt ~ 2 , where t is 
the lifetime of the FR II ( iKaiser & Besll2007l) . Since t oc p 1 / 3 , 
we have that p oc p 1 / 3 . The enthalpy is proportional to pV, 
but the volume of the lobe is an observable and is therefore 
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FIG. 20. — The estimated ICM heating rates due to FR lis (filled circles) and the estimated ICM cooling rates (solid line). The heating rates represent the 
enthalpy of each FR II in our sample divided by the age of the FR II. Since the enthalpy is calculated from the pressure inside the lobe rather than the external 
pressure, these points are upper limits on t he average heat i ng rat e. The cooling rate (solid line) is calculated by combining the scaling relatio n between X- 
ray lum inosity and velocity dispersion from Ortiz-Gil et al. (2004) with the scaling relation between velocity dispersion and cluster richness from Be cker et al.l 
1 2007). (A similar figure by Best et al. 2007 included several other scaling relations between X-ray luminosity and velocity dispersion; since these are all very 
similar to the Oritz-Gil relation, they are omitted here for clarity.) The jet powers (open circles) are higher than the required cooling rate, but only ~i% of this 
power generates the enthalpy of the lobe. Approxi mate ly 10% of the jet power is needed to balance cooling. (For comparison, the open circles represent the jet 
powers of the FR lis in our sample, estimated in 33.3| . The lower dotted line represents the median heating of FR lis in our sample assuming an uncertainty 
factor of / = 1 (i.e. assuming that our estimate of the ICM heating is exactly equal to the true ICM heating). The upper dotted line represents the median ICM 
heating assuming an uncertainty factor of / = 15. This is approximately the value of / required for FR II heating to balance the ICM cooling. 



independent of our assumptions about the ICM density. Thus, 
Hocp 1/3 . (13) 
If the density profile does not have a power-law slope of /3 = 
2, there will be an additional dependence of the enthalpy on 
the jet power, but if the slope is nearly 2, the deviation from 
the above relationship will be small. We therefore expect the 
dependence of the enthalpy on the ICM density assumption 
to be weak. If the typical ICM density is 5.9 x 10" 24 kg m" 3 
at 260 kpc (the largest density we use in this paper), rather 
than our assumed value of 7.2 x 10~ 26 kg m" 3 at 391 kpc, the 

enthalpy increases by only a factor of3. 

With a similar enthalpy analysis, iBest et ail (120071) found 
that radio-mode feedback has particular difficulty quenching 



cooling flows in very rich clusters. They estimated the aver- 
age ICM heating due to the AGN with a power-la w relation- 
ship b e tween radio lumi nosity and jet power from lBest et al.l 
(2006). Bes t et all J 2007 ) compared this heating rate to the ex- 
pected cooling rate estimated from several generic scaling re- 
lations between X-ray luminosity and velocity dispersion and 
found that the estimated heating rate fell below the cooling 
rate derived from the X-ray luminosity for clusters with a ve- 
locity dispersion above a v > 300 km s" 1 . All of the MaxBCG 
cl usters are expecte d to be above this velocity dispersion. 

IBest et aD (120071) identified three major sources of uncer- 
tainty in their analysis: (1) the conversion from radio lumi- 
nosity to jet power, (2) the estimate of the ICM heating from 
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the jet power, and (3) the fraction of ICM heating that takes 
place within the cooling radius. They parametrized these un- 
certainties by the ad hoc uncertainty factor, /, in EquationfTOl 
An uncertainty factor of / = 1 implies that the IC M heating 
within the cooling radius is exactly the amount that Be st et all 
(2007) infer from the radio luminosities. Larger values of / 
correspond to more ICM heating within the cooling radius 
than predicted, and smal l er valu es correspond to less heat- 
ing. Although iBest et all (120071) argue that / is likely close 
to unity, they note that if / is somewhat larger than 10 then 
the ICM heating due to radio-mode feedback would be suf- 
ficient to quench cooling flows in nearly all galaxy clusters. 
A better constraint on / could therefore determine whether 
radio-mode feedback is a viable mechanism to quench cool- 
ing flows in rich clusters. Because our approach takes ob- 
served FR II lengths and radio luminosities and uses an FR 
II evolution model to constrain the jet power distribution of 
our sample, we have effectively removed the first source of 
uncertainty fro m /. 

As noted by IBest et ail (120071) . there is additional uncer- 
tainty in the fraction of the heating that takes place within the 
cooling radius. We have assumed that all of it does, but this 
fraction could be much less than 1 . If this fraction is small, 
however, it becomes yet more difficult for radio-mode feed- 
back to quench cooling flows in rich clusters. 

Our res u lts sh own in Figure [20] are similar to those of 
IBest et ail (120071) — radio-mode feedback can only quench 
cooling flows in rich clusters (A^oo > 10) if the uncertainty 
factor satisfies / > 15. Meeting this criterion requires that 
the cores of galaxy clusters be much denser than expected, 
that FR II lobes be nearly spherical, and that the ICM be only 
slightly underpressured relative to the lobe. This strongly con- 
strains the environments in which radio-mode feedback is an 
effective mechanism to quench cooling flows. 

Even if / > 15, our duty cycle calculation suggests that 
FR lis can be an effective ICM heating mechanism for only 
a small frac tion of the cluster's lifetime. As described in de- 
tail in EU only -2% of the BCGs in the MaxBCG catalog 
have FR lis after the correction for selection effects. If every 
BCG in the MaxBCG catalog is as likely to host an FR II as 
any other (a reasonable assumption given the similarity be- 
tween the vast majority of MaxBCG clusters in redshift and 
richness), then FR lis are powered for ^2% of the cluster's 
lifetime. If the uncertainty factor / were greater than ~15, 
FR lis would be able to quench cooling flows during this ^2% 
of the cluster's lifetime, but would be unable to quench cool- 
ing flows at any other time. For FR lis to heat the ICM suf- 
ficiently to quench cooling flows for a substantial fraction of 
the cluster's lifetime would require the uncertainty factor to be 
yet an order of magnitude larger. Since we have largely elim- 
inated the source of uncertainty (1) of the conversion from 
radio luminosity to jet power, and the source of uncertainty 
(3) of the fraction of ICM heating that takes place within the 
cooling radius only acts to make / smaller, the only way to 
increase / is in the source of uncertainty (2) of the estimate 
of the ICM heating from the jet power. One way to increase 
/ is if there were heating due to shocks near the cluster cen- 
ter, as our calculations have only included ICM heating due 
to adiabatic expansion. 

Another possible solution that addresses both the appar- 
ently small value of / and the low duty cycle is if the FR II 
continues to affect the ICM after the radio source is no longer 
apparent. Hydrodynamical simulations by ? indicate that 
FR lis can substantially effect the cooling properties of the 



ICM for significant period of time (up to an order of magni- 
tude longer than the lifetime of the jet) after the jet has turned 
off by removing gas that otherwise would have cooled on 
longer timescales. Given a longer period of time in which to 
couple to the external gas, FR lis could be more efficient than 
their lifetime alone suggests. Our results therefore imply that 
if FR lis are major sources of heating in galaxy clusters, they 
contribute the vast majority of this heating after their demise. 

The range of jet powers of these FR lis provides some po- 
tential constraints on the jet powers in the additional 20% of 
the BCGs in the MaxBCG catalog that host other cluster ra- 
dio sources. The two most likely scenarios for why these 
other sources do not exhibit an FR II morpholo gy are that they 
fall b elow the critical jet power criterion (e.g. jKaiser & Best! 
120071) and that local variations in the ICM have di srupted their 
jets. With the notable exception of Cygnus A (Carilli et al. 
1994), nearly all of the early evidence for the impact of ra- 
dio sourcejijmjlte ICft£c as Perseus 
A (lBoehringeretal.1119931: iFabian et alJl2000l) and Hydra A 
dMcNamara et al.l[2000h . If the vast majority of cluster radio 
sources are not FR lis because they are insufficiently power- 
ful, then their jet powers must fall below the FR lis in our sam- 
ple, or < 5 x 10 36 W. This is at the low end of the estimates 
of the power requir ed to inflate cavitie s seen in X-ray gas of 
~ 10 35 " 39 W (e.g.. iBirzanetal .1120081) . although the clusters 
studied with X-rays may not be representative of all cluster 
radio galaxies. 

If the jet powers for these other radio sources are below 
the typical FR II jet powers by only one order of magnitude, 
or they are jets that exceed the critical value but have been 
disrupted, then they will be sufficiently powerful to counter- 
act cooling in their clusters. The first scenario would require 
the coupling between jet power and the ICM to be closer to 
100%. The second scenario could be constrained with de- 
tailed comparisons between the cavity powers and lobe en- 
thalpies of FR Is and lis. In addition, while our data suggest 
that the cavity power may be substantially less than the jet 
power for FR lis, a different relation between cavity and jet 
power may hold for FR Is. If they are more equal, then many 
FR Is may be too weak to create FR lis. Better estimates of 
the environmental conditions and intrinsic properties that pro- 
duce the FR I/II dichotomy would also provide valuable new 
constraints on the properties of the radio sources that likely 
ultimately produce feedback on the ICM. 

7. SUMMARY 

We estimate the typical lifetime, jet powers, and ages for 
a new sample of 151 FR lis identified in data from the 
FIRS T radio survey for clusters in the MaxBCG catalog of 
iKoester et all d2007al) . We use radio luminosities, projected 
length measurements, and the KDA FR II model to show that 
the typical FR II lifetime is 1.9 x 10 8 yr. We find that the 
distribution of FR II jet power in BCGs is best described by 
a log-normal distribution rather than a power-law distribution 
as has been previously assumed. We furthermore find that the 
jet power estimated for FR lis in our sample is lower than tha t 
of other samples (e.g.. lKaiser & Bestll2007t iBird et al.ll2008l) . 
The major uncertainties in the jet power and lifetime mea- 
surements are the density of the ICM and the axial ratio of the 
jets. Reasonable ranges in ICM density and axial ratio affect 
the lifetime estimate at the factor of two level. 

We examine the relationship between the properties of 
FR lis and the properties of BCGs and the clusters. We find 
that while BCG luminosity is correlated with FR II jet power, 
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other cluster properties such as BCG velocity dispersion and 
cluster richness are not strongly correlated with either the jet 
power or the FR II lifetime. These results suggest that the in- 
trinsic properties of FR lis are not highly dependent on their 
larger-scale environment. 

We evaluate the heating of the ICM from FR lis in galaxy 
clusters assuming adiabatic expansion of overpressured radio 
lobes. The time-averaged enthalpy of the lobes implies that 
ICM heating due to FR lis is smaller than ICM cooling by 
over an order of magnitude. Although the jet provides the lobe 
with a sufficient amount of energy to quench cooling flows, 
FR lis cannot be major contributors to ICM heating without 
a viable mechanism by which this energy can be distributed 
throughout the ICM. This stands in stark contrast to the total 
jet power, which exceeds the cooling rate by approximately 
an order of magnitude. If the jet power could couple with the 
ICM with ~10% efficiency, this would suffice to counteract 
cooling, albeit for only the ^2% of BCGs that host an FR II. 
The larger fraction of other radio sources in BCGs may be 
lower-power jets that are insufficiently overpressured to form 
an FR II morphology. These sources would require nearly 
100% coupling efficiency to the ICM and jet powers only an 
order of magnitude below our sample to balance cooling. 
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